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Abstract 

We present an asymptotic analysis of the effects of rapid rotation on the ground state properties 
of a superfiuid confined in a two-dimensional trap. The trapping potential is assumed to be radial 
and homogeneous of degree larger than two in addition to a quadratic term. Three critical rotational 
velocities are identified, marking respectively the first appearance of vortices, the creation of a 'hole' 
of low density within a vortex lattice, and the emergence of a giant vortex state free of vortices in 
the bulk. These phenomena have previously been established rigorously for a 'flat' trap with fixed 
boundary but the 'soft' traps considered in the present paper exhibit some significant differences, in 
particular the giant vortex regime, that necessitate a new approach. These differences concern both 
the shape of the bulk profile and the size of vortices relative to the width of the annulus where the 
bulk of the superfiuid resides. Close to the giant vortex transition the profile is of Thomas-Fermi 
type in 'flat' traps, whereas it is gaussian for soft traps, and the 'last' vortices to survive in the bulk 
before the giant vortex transition are small relative to the width of the annulus in the former case 
but of comparable size in the latter. 

MSC: 35Q55,47J30,76M23. PACS: 03.75.Hh, 47.32.-y, 47.37.+q. 
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1 Introduction and Main Results 

Since their first experimental realization in 1995, atomic Bose-Einstein condensates (BECs) have become 
a subject of tremendous interest, stimulating intense theoretical activity. The BE-condensed alkali vapors 
that are nowadays produced in many laboratories offer a spectacular level of tunability for several ex- 
perimental parameters, making them a favorite testing ground for many intriguing quantum phenomena. 
Among these is superfluidity, i.e., the occurrence of frictionlcss flow, previously observed in liquid helium. 
Atomic BECs offer a valuable alternative to the latter system, since experiments with dilute gases of 
ultracold atoms allow to test theories of superfluidity in much more detail. 

From a theoretical point of view an appealing aspect of the field is its sound mathematical foundation. 
The most commonly used model for the description of BECs, the so-called Gross-Pitaevskii (GP) theory, 
is now supported both by extensive comparisons with experiments [A"IIFel| and a solid mathematical basis 
[LSSYj . Indeed, substantial advances on the connection between GP theory and many-body quantum 
physics have been made in recent years, leading to a rigorous derivation of both the stationary |LSYli 
iLlSl IBCPY| and dynamic aspects jESYll IESY2I [P] of the theory. 

One of the most striking features of superfluids is their response to rotation of the confining trap. 
Typical experiments (see |Fel) for further references) start by cooling a Bose gas in a magneto-optical 
trap below the critical temperature for Bose-Einstcin condensation. If the trap is set in rotational motion 
the gas stays to begin with at rest in the inertial frame, but if the rotation speed exceeds some critical 
value (called f2 Cl below), vortices are nucleated and remain stable over a long time span. This behavior 
is in strong contrast with that of a classical fluid, which in the stationary state rotates like a rigid body 
and thus remains at rest in the rotating frame. 

In a rotationally symmetric trap, the main two experimentally tunable parameters are the rotation 
speed Q m t and the strength of interparticle interactions, written as 1/e 2 with e > in the sequel. The 
ground state of the system strongly depends on the relation between these two parameters. 

In this paper we study the ground state of a rotating Bose gas in the framework of the two-dimensional 
GP theorjjj. We consider the minimization, under a unit mass constraint, of the GP functional 

£ p G h P ys [*]= / dr - lAmt )V\ 2 + {V(r) - ^ t r*)\^ + , (1.1) 

Jr 2 I £ J 

with 

A rot := r2 rot A r = fl mt re#, (1.2) 

1 Such a description is justified if the trap almost confines the gas on a plain orthogonal to the axis of rotation or, on the 
contrary, the trap is very elongated along the axis, in which case the behavior can be expected to be essentially independent 
of the coordinate in that direction. 

2 The subscript 'phys' stands for physical, i.e., this is the functional in the original physical variables, as opposed to the 
rescalcd functionals introduced in Sect. 1.1 below. Units have been chosen such that h as well as the mass are equal to 1. 
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where f2 ro t is the rotational velocity, r = |r| with r <G R 2 is the distance from the rotation axis, and 
stands for the unit vector in the direction transverse to r. The confinement is provided by a potential of 
the form 



with k > 0. We restrict to the case 2 < s < oo (anharmonic, 'soft' potentials) and shall mainly study 
the Thomas-Fermi (TF) limit e — > 0. The case s = 2 is special in many respects because the potential 
— f2 2 ot r 2 due to the centrifugal force is also a quadratic function of r. As a consequence an upper bound 
is set to the allowed values of the rotational speed and different physics is expected in the regime where 
the centrifugal force nearly compensates the trapping force (see [2] and references therein, [LSY2] and 



We have studied before jCDYll ICRYl ICPRY11 |R2] the case of a flat trap with the unit disc as 
boundary. It corresponds formally to the potential 



which amounts to restricting the integration domain to the unit disc. This model can formally be regarded 
as the s —> +oo limit of the theory with the trap (|1.3[) but we emphasize that this statement has to be 
taken with care. Indeed, as revealed by the analysis of the present paper (sec also CPRY2 , Sect. Ill) 
the limit s — > +oo cannot be interchanged with the TF limit e — ¥ 0, which is the basis of most rigorous 
treatments of the GP theory, including the present one. 

Provided the confining potential is stronger than quadratic, i.e., s > 2, one can distinguish three 
critical speeds at which major phase transitions occur in the BEC (references in the literature concerning 
these phenomena include [UDJ |FBj EM E3 iKTUl iKBl iKFp : 

• When f2 rot < f^d the condensate is vortex-free. Vortices start to appear when the rotation speed 
exceeds tlr, and form regular patterns. An adaptation of the methods of |AAB1 IIM11 IIM21 lAJRj 
shows thatj f2 Cl tx e 4 /( s+2 )| loge|. This is valid for all traps, including the harmonic trap s = 2 and 
the flat trap, since the centrifugal force is negligible in this regime. 

• When Sl Cl <C S! rot <C Q C2 vortices are densely packed and uniformly distributed in the condensate. 
It is observed both numerically and experimentally that they arrange themselves into triangular 
arrays (Abrikosov lattices with hexagonal unit cells). It remains an open problem to provide a 
rigorous proof of this fact but a heuristic argument based on an electrostatic analogy and Newton's 
Theorem can be found in jCY] and [SS3] contains a rigorous proof for a simplified model. The 
distribution of vorticity has been shown to be uniform in this regime for the flat trap case both 
with Neumann [CYj and Dirichlet boundary conditions jCPRYlj . In this paper we adapt this result 
to the case of soft potentials. 

• When O ro t reaches the second critical speed il C2 , the centrifugal force dips a hole with strongly 
depleted matter density around the center of the trap. Vortices are however still uniformly dis- 
tributed in the annular bulk where the mass of the condensate resides. In the flat trap case the 
second critical speed is of order £ _1 [CDY11 ICRYi ICPRYlj . As shown in |CDY2) and discussed 
further in Section ITT21 below, this order is changed in the case of a soft potential (|1.3|) to e~~ , 
provided s > 2. Notice that in the limit s — > oo we recover the result for the flat trap case, despite 
the subtlety of the exchange of limits (see below). Here the behavior of the gas differs markedly 
from that in a purely harmonic trap where centrifugal effects start to compensate the trapping 
potential in the regime f2 ro t = 0(1). In this case the state of the condensate is well approximated 

3 Note that in the cited papers the GP functional is rescaled from the outset, as in Eq. 111.1211 below, and in the rescaled 
variables the critical velocity is uj c i oc | loge|. 



V(r) := kr s + 




(1.3) 




(1.4) 
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by a wave function in the lowest Landau level of the magnetic Hamiltonian corresponding to the 
first term in ([TTTjl jABNl lA"Bl ILSY2] and the GP theory ultimately loses its validity [LeS] . 

• When f2 rot is further increased above a third critical speed f2 C3 , a new phase transition occurs where 
vortices are expelled from the bulk of the condensate. All the vorticity resides in the central hole 
created by the centrifugal force, resulting in a giant vortex phase. One of the main results of this 
paper is that, for a potential of the form (|1.3|) with s > 2, Q C3 is of order e~ =+2 . This converges 
to e~ A when s — > +00, which is significantly larger than the value of fl C3 in the flat trap case. The 
latter is, to leading order in e, given by (2/(37r)) (e 2 | logel)^ 1 |CRY[ IR2j . The limits e — ► and 
s — > +00 can thus not be interchanged in this regime. As we explain further below, this is due to 
significantly different physical properties of the condensate in the two kinds of traps close to the 
transition regime. 

• When il Iot ^> Q C3 , the BEC is in a giant vortex phase and the energy is given to an excellent 
approximation by a trial function with a single multiple vortex at the origin (of degree f2 ro t to 
leading approximation), accounting for the macroscopic circulation. Nevertheless we prove that 
rotational symmetry is still broken in any true ground state, namely that no ground state can be 
an eigenfunction of angular momentum. 

The TF approximation in the regime f2 ci <C £l I0 t *C fi C3 was studied in |CDY2j where a preliminary 
upper bound for the vortex contribution was given. Precise upper and lower bounds are provided in 
Theorem 1 . 1 below. The case of large rotation speeds with fixed coupling constant e~ 2 has been considered 
before jRlj for the special case of a quartic potential (i.e., s = 4). The present contribution improves 
the results of that paper in several directions. In particular, although the appearance of a giant vortex 
state has been proved in |Rlj . no rigorous estimate of the third critical speed could be provided. Such 
an estimate is given in Theorem 1.3 below. 

In the remaining subsections of this introduction we define our setting more precisely and state our 
main results. We start by discussing the scalings we shall use to study the functional (jl.ll) . We then 
state our results about the regimes f2 Cl <C f2 ro t <C ^c 3 and f2 rot cx Q C3 in two different subsections and 
finally present our symmetry breaking theorem. 

For a concise discussion of our results and methods the reader is referred to the companion paper 



1.1 Scalings of the GP Functional 

In this paper we focus on the analysis of the minimization of the GP energy functional under the 
mass constraint 



The trapping is given by the potential (|1.3[) with £l osc < f2 ro t, and < e -C 1, i.e., we study the TF limit 
of the model. The power s in (|1.3[) characterizes the homogeneous trap and we assume that 2 < ,s < 00. 
We also introduce the parameter 



Since we are interested in exploring the rotation regime f2 ro t — > 00 but want also to keep track of the 
effect of the quadratic term in the expression above we shall assume that 



|OPRY2j . 





(1.5) 




(1.6) 



«eff = 7^< 



'rot 



(1.7) 
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for some given < 7 < 1. 

Since the profile of a minimizer of expands if f2 ro t — > 00 and/or e — > it is natural to rescale 
all lengths as well as the rotational velocity. In fact it is convenient to employ two different scalings in 
different asymptotic regimes: If 

n io t<e-^, (1.8) 
we define x € M 2 , x = |x|, ip(x), uj and A u by writing 

r = R E x, *(r) = fl-V(x), fl rot = R- 2 u, A w = wie fl , (1.9) 

where 

i? e := (fc e 2 )~^2 =C (e'^Y (1-10) 
With these scalings the functional (|l.ljl can be rewritten as 

£ p G 4[*]=i?- 2 ^ p ^], (1.11) 

with 

£S P m-.= f dx{i|(V-iA u )^| 2 +e- 2 [x 8 M 2 -h^ 2 ^H 2 + W 4 l}- (1-12) 
If on the other hand 

aot>e"Si, (1.13) 

a natural scaling parameter is given by the position R m of the unique minimum point of the effective 
potential (jl.6p . which is explicitly given by 

/"rot \ I 



Writing now 

r = i? m x, tf(r) = fl"V(x), a-ot = i? m 2 ^ A n = fia«*, (1.15) 
the GP energy functional (|1.1[) becomes 

fp G h p y s[*] = i? m 2 [eSTM + (i - 1) 7« 2 ] , (1.16) 

with 

£§ P bP}~ [ dx{i|(V-zAo)^| 2 +7^(x)|^| 2 + £ - 2 Ki 4 }, (1.17) 
where the effective potential is given by 

t s — 1 t 2 — 1 

W(*) := ^ - (1-18) 

Note that for convenience we have extracted a negative additive constant in (|1.11|) exploiting the L 2 - 
normalization of both ^ and ip, in order to obtain a positive trapping potential W > with 

inf W(x) = W(l) = 0. (1.19) 

Note also that 

O rot = (k 2 )^w = {sk/-f)^2ft^ (1.20) 
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and the two rescaled velocities, w and fi, are related by 

UJ = (s/7)^T5£~^T3ftSf . (1.21) 

For s — > oo, i.e., a fiat trap, u = Q = fi ro t- 

To find the ground state energy and wave function of the system, we minimize the functionals (|1.12D 
or (|1.17p on the domain 

^ GP := {V € H^R 2 ) : |jV|i 2 = 1} . (1.22) 

We introduce the notation 

E^ P := inf^ P M, ^ P := , i£ £n P [# (1-23) 

for the ground state energies. Any corresponding minimizcr will be denoted by tp^ or respectively 
The GP minimizcr is in general not unique because vortices can break the rotational symmetry (see 
Section ll.4[) but any minimizer satisfies in K 2 a variational equation (GP equation) which is in the case 
of (fl~T2jl 

- iA^ p - uL z ^ + a-* [r° + 1(1 - ^eW + 2 |^T] ^ p = ^ P ^ P , (1.24) 
whereas for (| 1 . 1 T[) it becomes 

- ±A^ P - n Lz ^ p + ifiV^f + l&W{x)^ + 2e-z \^ p \ 2 i>g p = MH P V'g P . (1-25) 

Here L z = —idtf stands for the z-component of the angular momentum and the chemical potentials 
and //q P are fixed by imposing the L 2 -normalization condition on tp^ F an d "0o P respectively and are 
given by 

MS P = ^ GP +W dx|^ p | 4 , /4f^ P +W dxUg P | 4 . (1.26) 



The functional ()1.12j) is used to analyze the behavior of the supcrfluid when fl Cl f2 rot <C fi C2 1 which 
corresponds to |loge| <C u> <C £ _1 whereas the scaling leading to (|1.17[) is more suited to the regime 
f2 r ot 3> ri C2 ! corresponding to e _1 <C The transition to a giant vortex phase occurs when f2 oc e~ 4 . 

For the sake of simplicity we will often omit the suffix ui or in t/i gp \ ■> , etc., when it is 
clear the regime under discussion and/or the statement applies to all the regimes taken into account in 
the paper, i.e., | loge| <c oj, il < e~ 4 . 

1.2 Between Q Cl and f2 C3 : Vortex Lattice and Emergence of a Hole 

As we are going to see the leading order contribution to the GP energy E is determined in this regime 
by the density p GP := |?/> GP | 2 , not the phase of the order parameter. Correspondingly we introduce a 
TF-likc functional acting only on a density p: If 

| log e | < w < e-\ (1.27) 

we set 

<TM = ^ / dx [(X s +p)p- ±jeWx 2 p] , (1.28) 

JR 2 

whith p belonging to the domain 

:= |p G L 2 (M 2 ), \\p\\i = 1 : p > 0,x s p e L^M 2 )} . (1.29) 
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The minimization of £ TF is straightforward (see Section 12. ip and both the energy and minimizer have 
explicit expressions: The latter is the compactly supported function 

pY(x) = \ -* S + freW] + . (1-30) 

where /i TF is determined by the L 1 -normalization and [ ■ ]+ stands for the positive part. We also set 

i£ F := M8™[p}=£r[p T f]- (1-31) 
pes> if 

On the other hand i£| 

£- 1 <0<£- 4 , (1.32) 
the scaling defined in (|1.15[) becomes more useful and it is convenient to define the TF functional as 



£^[p} = f dx[e- 2 p 2 + 7 n 2 W(x)p] 
Jr 2 



(1.33) 



with p € ^ TF . The TF energy and minimizer, defined analogously to (|1.31[) . are now denoted by Eq F 
and pj^ respectively. The explicit expression for the latter is 

p^(x) = ±[e 2 ^ -je 2 n 2 W(x)} + . (1.34) 

As we are going to see (see Section |2~T1 for further details) there are critical values oj c , S7 C = 0(e^ 1 ) of 
the angular velocity satisfying 

Sl c = R 2 m R- 2 LJ c , (1.35) 

such that if u) > ui c (rcsp. S7 > Sl c ) then the TF density p^ F (resp. pj^) vanishes at the origin. More 
precisely, if SI > Sl c , there exists a disc of radius x- m > so that 

pJi \ x ) = 0j f° r an y — x — x in- (1.36) 

We also point out that by the rescalings (jl.9|) and (|1.15|) . lu < lj c holds if and only if SI < Sl c , so that the 
conditions | loge| <C uj < lj c and Sl c < SI e -4 cover all the possible asymptotic regimes between Sl Cl 
and S7 C3 . 

We can now formulate the result about the energy asymptotics of E GP : 

Theorem 1.1 (GP ground state energy). 

// | loge| uj < oj c as e — >• 0, then 

ES P = EZ ¥ + \u |log (A.) I (1 + (1)), (1.37) 

while, i/Sl c < SI <C e~ A , 

E£ P = EE F + ±fi |log ( £ 4 S1) | (1 + o(l)). (1.38) 

As we have anticipated a transition takes place when the angular velocity gets larger than some 
critical value of order 0(e~ 1 ): A direct consequence of the asymptotics (|1.37|) and (|1.38[) (see, e.g., 
[CPRY11 Proposition 3.1]) is that 

lll^T-P^lla^oWllp^L, (1-39) 

i.e., the density |i/j gp | 2 is close in L 2 (M 2 ) to the TF density p TF , but when SI gets larger than Q c the 
latter develops a hole and a macroscopic region where ?/> GP is exponentially small in e appears close to 
the center of the trap. Next proposition shows that the threshold of such a transition is precisely given 
by Sl c , i.e., 

S1 C2 = (**/ 7 )*n^, Sl c = 2^|^ (!) " e \ (1.40) 



4 Notice that by the scaling relation 111.211 1. Q = C(£ _1 ) if and only if w = 0{e~ 1 ). 
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Proposition 1.1 (Emergence of a hole in i/j GP )- 

If, as e — > 0, > r2o£ — 1 with Qq > efl c there exists < Xhoic < 1 such that 

|V GP (x)| =<D(s°°), (1.41) 

for any x such that < x < Xhoic • 

We conclude this section by stating a result about the vortex distribution and showing that throughout 
the parameter region |loge| < lo < e _1 (resp. e" 1 < 17 <C £~ 4 ), the vorticity of ip^j p (rcsp. -0n P ) is 
uniformly distributed on a length scale S7 -1 / 2 inside supp (/oJ F ) (rcsp. supp (/Qq F )) (see also the Remark 
after Theorem II ,2|) . For technical reasons, however, we are not able to prove the result in the whole 
support of p TF , but we have to restrict the analysis to the slightly smaller region 

ftbulk := {x G R 2 : x< < x < x> } , (1.42) 

where for some positive parameter /3 := /3(e, f2) such that f3 — >• as e — s- (see (|2 . T6|) for its precise 
definition) 

^x in + /3(efi) ' , if f2 c < O < £ 4 , 

Here a;i n and x ou t are the inner and outer radii of the support of the TF density. The definitions (|1.43l) 
imply that 

|su PP ( P TF ) \ TZ hulk \ = O (j8(l + (en) 2 ^)- 1 ^ < |supp (p TF ) I , (1.44) 



and 



\^ GP \\r^, ^=l-o(l), (1-45) 



li 2 (K bu ik) 

i.e., 72-buik contains the bulk of the mass of ip . 

In the following theorem, we describe the vorticity distribution in the bulk of the condensate for 
rotation speeds between the first and third critical speeds. As it is usual in the Ginzburg-Landau (GL) 
theory, this is done by investigating the asymptotic properties of a certain vorticity measure whose 
definition is part of the statement of the theorem. Note that we do not tackle the question of the precise 
location of the vortices, which should be on the sites of a hexagonal lattice. Proving this fact rigorously 
remains a major open problem in both GP and GL theories, see |SS3] for the most recent results in this 
direction. 

In the statement of the theorem we use the notation £>(x, p) for a closed disc ('ball') with center at x 
and radius p. 

Theorem 1.2 (Uniform distribution of vorticity). 

If | logej <C u < £ _1 or e^ 1 < -C £~ 4 as e — > 0, there exists a finite family of disjoint balls {£>;} := 
{B(xj, Qi)} C T^buik, i = 1, ■ ■ ■ ,N, such that 

1. Qi < o K 1 / 2 + n- 1 / 2 ) , £ f?? < (1 + (eO) 2 / 3 )" 1 ; 

2. |V> GP | > on dB t , i = l,...,N. 

Moreover, setting d{ := deg{"0 GP , dBi] and defining the vorticity measure as 

N 

v := 27r^d 4 (5(x-x l ), (1.46) 

i=l 

then, for any set S C 7£buik such that \S\ 3> log(£ 4 i7)| 2 as e — > 0, 

^1 -> 1. (1.47) 

Q\S\ e^O { ' 
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Remark 1.1 {Vortex balls) 

By a closer inspection of the proof of the above Theorem ll.21 one can realize that there is more information 
contained in the family of balls {£>;}: Each ball might contain a very large number of vortices but almost 
all vortices are enclosed by those balls. More precisely, there exists a 'good region' in the condensate 
where all the possible zeros of ip are contained in U£>; and outside this set ip GP is vortex free. There 
is also a complementary 'bad region' where we are not able to prove such a thing, but a key point of the 
proof of Theorem 11.21 is that the good region covers most of the condensate, i.e., the ratio between the 
areas of the bad and good regions is o(l) as e — s- 0. 



1.3 The Regime f2 oc e 4 : Transition to a Giant Vortex State 

In this section we consider f2 given as 

n = J < L48 ) 

where S7o is a fixed constant while e — > 0. 

A 'giant vortex state', that turns out to be an excellent approximation to the ground state for large 
Clo, is obtained by minimizing the GP functional over functions of the forrrj^ /(x)e 4 L f2 J 1? with / real. The 
giant vortex energy functional is given by 

:=3P[/e iLnJ *]= / dx{i|V/| 2 + 2 ^)/ 2 +e- 2 / 4 }, (1-49) 

where 

U(x) := \\n\ 2 n- 2 x- 2 + \x 2 + ± 7 {x s - 1) - ± 7 (x 2 - 1) - [n\Q-\ (1-50) 

By standard arguments (see Proposition 13. 3p one can show that there exists a unique minimizcr g gv of 
(|1.49[) in the domain 

^ := {f e H^R 2 ) : f = f* > 0, X s f G L 1 ^ 2 ), ||/|| 2 = 1} . 
In addition it is easily seen that <? gv is radial and does not vanish except at the origin. We also set 

E sv := inf £ gv [/l. (1.51) 

The leading order of the GP energy and density will be given with very good accuracy by the minimization 
of the simpler giant vortex energy functional. It is thus of importance to understand the properties of the 
latter. When ft e -4 , its minimization leads to a profile that is to leading order of TF-type as discussed 
in the preceding subsection. A salient feature of the regime (|1.48p is the fact that, when Slo gets large, 
the giant vortex density g gv changes from a TF profile to a gaussian distribution centered on the circle of 
radius 1. This is most conveniently understood by using an auxiliary one-dimensional energy functional 
(recall that g gv is radial) 

£^[f] = 2n [ dx{\\f'\ 2 + \ti 2 a 2 {x-l) 2 f 2 +e- 2 f} = 

2nQ [ d y{i|/'| 2 + iaV/ 2 + Oo 1/2 / 4 }, (1-52) 

where 



a := v/4+ 2(s - 2). (1.53) 



' Here [ J stands for the integer part. 
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One obtains £ aux from (|1.49[) by expanding the potential around its minimum at x = 1 and neglecting 
terms beyond quadratic in the expansion. The variable transformation y = Q 1 ^ 2 (x — 1), f(y) = f^ 1 / 4 f(x) 
is then used to blow-up to the characteristic size of the domain where the density is confined. It is clear 
that all three terms in (|1.52[) are of the same order of magnitude but the importance of the last term 
diminishes with increasing D,q. The quadratic part of the energy is simply the energy functional of the 
harmonic oscillator, so without the last term the minimizer is the Gaussian 

9osc(y) = ^ 1/4 « 1/4 cxp{-im/ 2 }. (1.54) 

We expect and prove (see Section I3T21 that g gv is close (after rescaling) to this function when Hq is 
large. This behavior of the density profile is very different from that given by the TF theory discussed 
before. The fact that the transition to a gaussian profile occurs in the same parameter regime as the 
transition to a giant vortex phase is at the origin of several new features of the physics of the third phase 
transition in a soft potential as compared to the theory in a flat trap. This does not seem to have been 
noticed in the literature before: As far as we know, the profile is always assumed to be of TF type in 
previous papers dealing with the transition to a giant vortex state, with the exception of [Rlj . 

To begin with, it is more delicate to define the bulk of the condensate in the regime large. The 
gaussian mass distribution has much slower decay and longer tails than a TF-likc profile. As a consequence 
the mass is much less concentrated and it is necessary to consider an annulus of thickness much larger 
than the characteristic scale Q.^ 1 / 2 = fl 1 ^ 2 e 2 of the gaussian distribution: 

Auik :={x£l 2 : 1 - c| loge| 1/2 fr 1/2 < x < 1 + c| logel 1 / 2 ^- 1 / 2 } , c < yj^. (1.55) 

Our main result about the regime (|1.48[) is the following: 
Theorem 1.3 (Absence of vortices in the bulk). 

If SI is given by (|1.48[) . there exists a finite constant f2o such that for any f2o > n o minimizer ip GP 
has a zero inside .Abulk if e is sufficiently small. 
More precisely, for any x £ -4buik 7 

|^ GP (x)| =g gv (x)(l + 0(\loge\- a )) (1-56) 

for any a > 0. 

The estimate (|1.56[) justifies the notation .4buik: Indeed we have 

f dx |V^ GP | 2 = l-o(l), 

J -Abuik 

because such an estimate holds true for g gv , so .Abuik contains the bulk of the mass. Theorem 11.31 then 
characterizes the transition to a giant vortex phase where the bulk of the condensate resides in an annulus 
free of vortices but with a macroscopic circulation around the origin. Recalling the results of Theorem 1 1.2 1 
we deduce that the third phase transition occurs when f2 is of order e~ 4 . Indeed, for speeds much smaller 
than this threshold, vortices do occur in the annular bulk, whereas they are absent if SI is a sufficiently 
large multiple of e~ 4 . 

The core of the proof of Theorem 11.31 is a detailed analysis of the asymptotic behavior of the GP 
ground state energy, leading to the following estimate: 

Theorem 1.4 (Ground state energy asymptotics). 

If O is given by (|1.48[) with S7q > CIq as in Theorem \1.SX then as e — > 



E QT = £gv + 0(| log£ |9/2)_ 



(1.57) 
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We are also able to estimate the degree of the GP minimizer, ensuring that there is a macroscopic 
phase circulation around the central low-density hole. This is the content of the following theorem. 

Theorem 1.5 (Asymptotics for the degree). 

// fi is given by (|1.48p with fio > Oo as in Theorem \1.3\ and R is any radius satisfying 

R=l + O(0.- 1 / 2 ), (1.58) 

then as e — > 

deg (V GP , dB R ) =n + 0(fi o | loge| 9/4 ). (1.59) 

As already mentioned, the physics of the giant vortex phase transition in a homogeneous trap turns 
out to differ markedly from the corresponding regime in a flat trap that was studied in CRYl [CPRY1 1 |R2] . 
As a result, the technique of proof is different and does not rely on any construction of vortex balls. The 
remaining part of this subsection aims at elucidating the main steps in the proof. To bring out the salient 
points technical details are skipped and most of the following equations will be written in a simplified 
form without the remainder terms. We remark also that the method to be presented is simpler than that 
in |Rlj and leads to sharper results, although some ideas are common to both papers. 

As a starter let us note that, in the soft potentials (s < oo) we consider, the centrifugal force does not 
squeeze the condensate against a hard wall as in the flat potential (s = oo) considered in |CPRY1| although 
it still concentrates the scaled density in an annulus around x = 1. The soft potential, however, allows 
more spreading of the matter density than a trap with a hard wall, so the density is lower and vortices 
therefore energetically less costly. This explains the fact that the third critical speed is considerably 
larger in a homogeneous potential than could be expected from the analysis of the flat trap model. 

It had been recognized before |Rlj that for extreme rotation speeds the density profile would be 
gaussian, but when studying the giant vortex transition a key difficulty originates from the fact that the 
change from the TF profile to the gaussian distribution takes place in the same parameter regime as 
the one where vortices are expelled from the bulk of the condensate. In a sense two distinct transitions 
happen at the same time in the regime (| 1 .48[) and this is the main reason why the analysis we perform 
in this paper is rather different from |CRY1 ICPRYlj . 

Indeed, in CRY, CPRYl] the potential vortex cores were much smaller than the width of the annular 
bulk of the condensate and vortices disappeared for energetic reasons. To analyze this phenomenon 
precisely we had to adapt the method of vortex balls originally developed in the GL theory (see |SS2j 
and references therein) and used to evaluate il Cl [IMli |A"JR] . 

In a homogeneous trap the cores of potential vortices turn out to be comparable in size to the width 
of the annular bulk when 51 ~ Q Ca . It is thus not clear whether they disappear because of their energy 
cost or are simply expelled from the condensate because they are too large to fit in. This explains why 
we have not been able to give a precise estimate of the third critical speed: Potential vortices are too 
large to treat them using the GL technology. Their energy cost can not be evaluated as accurately as in 
[CRY! IR2] and we have to rely on coarser estimates. 

The starting point and the basic methodology is however reminiscent to that employed in |CRY[ 
ICPRYlj and in works dealing with the first critical speed in GP theory jAABl ITMTl HM2] . We first 
decouple the energy functional, an idea originating in |LMj and used repeatedly in the literature. Actually 
it is more convenient to restrict the integration domain before performing the energy decoupling. For 
technical reasons it is not sufficient to work on .Abuik (defined by Eq. (|1.55p ) and we consider instead a 
larger domain 

A 7h = {x£l 2 : 1 - i le e 2 < x < 1 + r) £ e 2 } (1.60) 

where we have some freedom in the choice of ?y £ 3> 1. It should in any case satisfy the condition 
?y £ » | log el 1 / 2 in order for A„ s to contain .Abuik but since we shall also need to consider A^ we even 
require n £ 3> |loge|. If r\ e is chosen large enough, both the true GP minimizer and the giant vortex 
density profile arc exponentially small in e outside A„ E and one can restrict the analysis to this annulus 
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at the price of negligible remainders. On the other hand, the proof of energy estimates and the absence 
of vortices requires that rj E is not too large and we shall eventually choose r\ e = | loge| 3 / 2 , cf. Eq. (|3.104j) . 

By defining a giant vortex functional reduced to A„ E and minimizing it, we obtain an energy and 
a profile g„ E supported in the annulus A Ve ■ We then write any test function ip as 

$(*)=g^{x)e i M*v(x), (1.61) 

which defines v on A^ ■ Since g riE is without zeros, the function v contains all possible zeros of ij) in the 
annulus. The variational equation for g 11e leads to 

£ GP ^}^E^+£ nE [v] (1.62) 

with 

£ Ve [v]= f dxgl{l\Vv\ 2 -B-(iv,Vv)+e- 2 gl(l-\v\ 2 ) 2 }, (1-63) 
where B = (tlx — [^Ja; -1 ) and (iv, Vf) = ^(uVu* — v*Vv). It is plausible that if we write 

V' GP (x) = .g, )e (x)e lLf2J ' 5 u(x), (1.64) 
then the function u defined on A, h should minimize £ n _ under the constraint 



/ dx<4M 2 = l. (1.65) 



This is not fully rigorous because some remainder terms in (11 . 62[) are neglected, most of which can be 
controlled using the exponential smallness of ip GP outside of A„ E ■ Nevertheless, studying the reduced 
problem of minimizing £ Ve under the constraint f| 1 . 65[) gives some valuable insight on the full proof we 
are to develop in Section [3] 

The appropriate trial function is given by v = 1, which corresponds to a pure giant vortex ansatz. 
This gives 

£ 1h [u] < 

and the gist of the proof is the obtainment of the corresponding lower bound, proving that the giant 
vortex trial function optimizes the energy up to small remainders. To this end one needs to control the 
only possibly negative contribution originating in the second term of (|1.63[) . This is most conveniently 
done using a potential function defined as in |AAB[ ICRYj : 

f(x)= f dt (fit- Lfijr 1 ) (1-66) 



This function satisfies V^i 7, = ff^B and vanishes on the inner boundary of A„ E . Using Stokes' formula 
on the second term of (|1.63p leads to 

£„.[«]=/ dx{i ffr 2 ,jVu| 2 + J F Ai + £- 2 4 e (l-| U | 2 ) 2 }-^(l + e 2 % ) f da(iu,d T u). (1.67) 

J Arte JdB l + .^ s 

The quantity fi is defined as 

H := curl(Mi, Vit) (1.68) 

and is often referred to as the vorticity measure because in many situations involving vortices one can 
prove that such a measure is close to measure of the form (|1.46l) : a sum of Dirac masses counting the 
vortices of u with multiplicity. It is not completely clear that fi deserves such a denomination in the 
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regime we consider: As already emphasized, the cores of potential vortice^l are too large to employ the 
jacobian estimate of |JS| that would put the comparison with a sum of Dirac masses on a rigorous ground. 
For illustration we shall employ the terminology nevertheless. It is justified by the fact that only the 
term involving \i in £ r)e [u] could favor the nucleation of vortices. Indeed the other bulk terms are positive 
and the boundary term should morally be negligible as we now explain. 

Since Q is very large it is plausible that one can replace it with its integer part at negligible cost. 
Also the density profile should be almost symmetric about the radius 1. If it were exactly symmetric, 
F(l + 7] s e 2 ) would be extremely smalfl and the boundary term in (|1.67p could be neglected, reducing the 
proof of the lower bound to the estimate of the second term in (|1.67|) . 

To deal with the term f , Fu we first remark that due to the Cauchy-Schwarz inequality we have 

\fj,\ < \Vu\ 2 (1.69) 

pointwise, and in the regime (|1.48[) g Ve and F are of the same order of magnitude. It is thus plausible 
that g^jVit| 2 can control F^i pointwise. This is in marked contrast to the situation in a flat trap where 
g Ve is much smaller than F, making the method of vortex balls indispensable to improve the precision of 
the estimates. 

For fio large enough, g, h is approximately gaussian, so that one can use the ansatz () 1 . 54[) to compute 
F and obtain 

\F{x)\<a- x g 2 Ve {x) (1.70) 

for any x <G -4buik and large enough. As a > 2 (see (|1.53|l ). combining (|1.69p with (|1.70[) we can 
conclude that the kinetic energy term in (|1.67[) dominates the term involving the vorticity measure, which 
leads to a lower bound confirming that the giant vortex ansatz is optimal (up to negligible remainders). 
As a by-product of the proof one can also obtain an estimate of the third term in (| 1 .63(1 that contradicts 
the presence of vortices in >4buik: The technique, involving an estimate of ||Vu|| Loc , originates in |BBHj 
and has already been used in the flat trap case |CRY[ ICPRYlj . 

What remains to be understood is the way in which one can justify the approximations made above. 
We have assumed that f2 could be replaced by its integer part. This is harmless and the corresponding 
remainder terms are ignored in what follows. More importantly we have assumed that 

• g r)c is a radial Gaussian centered on x = 1 for f2o large enough, which allows to prove (jl.70[) : 

• in particular g Vs is symmetric with respect to 1 so that the integration by parts in (|1.67[) produces 
only very small remainder terms. 

It is fairly easy to prove that g Ve is close to a Gaussian around the center of the bulk profile where the 
matter density is large enough. However, we need to justify the pointwise estimate (|1.70[) in a domain that 
is much larger than the bulk, including a region where the density g Ve is extremely small. To obtain an 
appropriate lower bound to , a subsolution is constructed for the corresponding variational equation. 
For the method to apply it is necessary to restrict to a smaller annulus: We are able to prove that 
is pointwise close to the optimal Gaussian only in A^j^- The region A Ve \ A^: is dealt withal using the 
exponential smallness of ijj GP m the low density hole. 

As for the boundary term in (jl.67[) one can prove that g Ve is to a very good approximation symmetric 
around x = 1. Indeed, as we have explained, the profile becomes gradually gaussian when Qq is increased, 
because the weight of the quartic component of the energy decreases. The reason we have separated the 

6 The size of vortex cores can be estimated by taking it as a variational parameter in a trial function with prescribed 
zeros and minimizing the sum of the first and third terms in 111.6311 . 

7 As t — t^ 1 is to a very good approximation antisymmetric with respect to 1, F(l + rj £ e 2 ) would reduce to the integral 
of an odd function on a domain symmetric with respect to the origin, which vanishes. 

8 Thcre is some freedom in the choice of A.^rr, the essential point being that we can prove the appropriate lower bound 
on g r]e only in a smaller domain than A\ le ■ 
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discussion of the symmetry from that of the approximation by a Gaussian is that g, h is very close to 
a symmetric function even for fio small where the gaussian approximation is not expected to be valid. 
In some sense g^ is much more symmetric than it is gaussian. To understand this, one can consider 
the energy functional (jl.52[) obtained by keeping only the quadratic term in the Taylor expansion of the 
potentials: After rescaling, the minimizer is symmetric for any £Iq and one can prove that it approximates 
g Ve much betteJHI than the Gaussian (|1.54[) . 

Making use of this property we show that F(l + e 2 r) £ ) is indeed very small. Then one could hope to 
conclude that the boundary term is negligible if a proper estimate of the phase circulation 



were available. This is however not easy to obtain because the integral is located in a region where g ne 
is very small, casting doubt on the possibility to prove any estimate of the term above using the kinetic 



To get around this point we modify our strategy as follows: Instead of the potential function (|1.66ll 
we use two different potentials F\, F% vanishing respectively on dBi_ e 2 rje and dBi +e 2~ . The first is 
used to calculate the part of the momentum term residing between the inner boundary of A Ve and the 
maximum point of g, u , while we use the second in the complementary region, leading to a formula of the 
form 



with the convention that F\ = for r > 1 and F% = for r < 1. In some sense we are using a 
discontinuous potential, which has the effect to produce a boundary term on dB±, i.e., in the middle of 
the bullf^l. An important point is that I-F1 + -F2I < \g\ c when flo is large, which means that the analysis 
of the bulk terms can be carried out exactly as sketched above. 

On the other hand the jump of the potential function -Fi(l) — -^2(1) is roughly speaking given by 
F(l + e 2 i] £ ) with F defined in (|1.66[) and is thus very small because of the symmetry of g Vc . The key 
point is that this information can be combined with an efficient estimate of the phase circulation of u on 
dB\. Indeed, the boundary term is now located where the matter density is large, which allows to control 
the phase circulation of u in terms of the kinetic term in the energy (see Lemma !3-5[) . This control also 
leads to the degree estimate of Theorem 11.51 once u has been proved to be vortex free. 

The trick of using a discontinuous potential function to make boundary terms tractable in the study 
of functionals such as (|1.62|) seems to be new. In situations considered before |AAB1 ICRY1 ICPRYlj , the 
TF-like profile allowed to employ other methods for estimating boundary terms. These fail in the context 
of the present paper due to the long tails of the approximately gaussian density profile. 

1.4 Rotational Symmetry Breaking 

The nucleation of vortices in a rotating superfluid trapped in a radial potential has an important connec- 
tion with the symmetry properties of the ground state. For example, when il rot is kept below the first 
critical speed, the ground state is expected to be unique and rotationally symmetric. This has recently 
been proved in [AJR] at least in the TF limit: The ground state for f2 rot < fl Cl is equal to the ground 
state without rotation, i.e., for f2 ro t = 0, which is a radial function. When there is only one vortex in 
the superfluid, it will sit at the center of the trap and the ground state, although not radial any longer, 
is still an eigenfunction of the angular momentum, i.e., it is of the form f(r)e with / a real-valued 
function and d £ Z. 

9 One can compare the remainders in 113.431 and 113. 6711 below. 

10 More precisely the jump is at the maximum point of the profile, which is close to 1 but not exactly equal to it. 




energy \ g\ 




£ V M- I te{yi\Vu\ 2 + (F 1 +F 2 )» + e- 2 gl(l-\u\ 2 ) 2 }-(F 1 (l)-F 2 (l)) [ da(iu,d T u), (1.71) 
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On the other hand the occurrence of more than one vortex in the GP minimizer 4> is automatically 
associated with a breaking of the rotational symmetry of the functional at the level of the ground state 
(see, e.g., jSeirl ICDY1] for a detailed discussion of this phenomenon). Therefore Theorem 11.21 guarantees 
that rotational symmetry is actually broken for any | loge| <C u), -C £~ 4 , since the bulk of the condensate 
is filled with a large number of vortices, namely of the order of the angular velocity lo or f2. 

However, the giant vortex transition described in Section 11.31 seems to suggest a restoration of the 
symmetry under rotations: For instance the proof of Theorem 11.41 shows that a giant vortex trial function 
with all the vorticity concentrated at the origin can reproduce the GP ground state energy up to sublcading 
order corrections. In addition the bulk of the condensate is vortex free (see Theorem 1 1.3[) so all vorticity 
must reside in a region where the density is low. 

In the following theorem we prove that the rotational symmetry is never restored even though the 
angular velocity crosses the giant vortex threshold. This suggests that there is a nontrivial distribution 
of vortices in the low density hole and it is an interesting open problem to locate them, or else find other 
reasons for the symmetry breaking. 

Theorem 1.6 (Rotational symmetry breaking). 

If Vt > £~ 4 , no minimizer tp GF of the GP functional is an eigenfunction of angular momentum for e 
sufficiently small. 

Note that the theorem is stated for f2 > e~ 4 , i.e., in the giant vortex regime. This is the regime where 
one could a priori hope that the ground state of the theory is exactly of the form f(r)e lM . 

The same result has already been proven for smaller angular velocities in |Seirl Theorem 2]: Assume 
for instance that | loge| < u> <C e _1 , then it is not difficult to realize that the total winding number of 
any minimizer must be O(to). Therefore if one picks some finite d > s/2 + 1, Theorem 2 in |Seirj proves 
the instability of any symmetric vortex of degree n = O(w): Condition (2.29) in |Seir| is indeed always 
satisfied for e sufficiently small, since the chemical potential is of order 0(e~ 2 ), whereas the r.h.s of (2.29) 
is 0{lu 2 ). 

1.5 Organization of the Paper 

The proofs of the results stated above are given in the following sections. The vortex lattice regime is 
discussed in Section 2, beginning in Sect. 2.1 with a comparison between the TF energy and density with 
those of an energy functional where the vector potential responsible for the appearance of vortices has 
been dropped. Section 2.2 contains the proof of upper and lower bounds to the vortex contributions to 
the energy. The uniform distribution of vorticity is proved in Section 2.3. 

Section 3 is concerned with the giant vortex transition. After some preliminary estimates on the 
ground state energy and minimizing wave functions in Sect. 3.1, we study in Sect. 3.2 the giant vortex 
density profile. The essential results here are the estimates on the deviation from a Gaussian (Proposition 
3.5) and from another profile symmetric about x = 1 (Proposition 3.7). The basic energy estimates are 
given in Sect. 3.3 and in Sect. 3.4 they are combined with gradient estimates to establish the giant vortex 
transition. 

Finally, breaking of rotational symmetry is proved in Section 4. 

2 The Regime | loge| « w, O « £ 4 
2.1 The TF Energies and Densities 

We start by summing up some important information about the minimization of the TF functionals £^ F 
and £^ F (Equations p.28[) and (|1.33l0 : We often refer to |CDY2] for a detailed analysis and omit here 
the proof details. 
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The unique TF minimizer of the functional f| 1 . 28[) for | loge| <C u> < e 1 is the density (|1.30|) . i.e., 

with fi^ F = E TF + £~ 2 \\Pu F \\2- The support of p^ F is always a compact set since s > 2, and there exists 
a critical value of the angular velocity uj c = 0(e~ 1 ) such that 

/ tf\ J {x G M 2 : x £ [0, x ut]i , ifw<w c , 
supp ( Puj = i } J (2.1) 



with Xin > 0. It is given by 



If a; = u) c , the outer radius is 



For generic w neither Xi n nor x out are explicitly given but can be obtained by solving the equations 
/°3 F ( a; in/out) — 0- Both radii are increasing functions of u, whereas both p^ F and x ou t — x- ln . i.e., the 
width of supp (p3 F )j decrease as uj increases. 

When uj becomes larger so that uj > e _1 (which is the same as f2 > we use the definition (jl.33l) 
for the TF functional. The TF minimizer is in this case 

pZ F (x):=l[e 2 t£ F - 7 eWW(x)} + . (2.3) 

The support of pJ F is always compact, but if O > e" 1 is sufficiently large, the density vanishes in a 
neighborhood of the origin. More precisely one can define an analogue, denoted by f2 c , of the critical 
value uj c as the smallest value of fi such that pJ F vanishes at the origin x = 0. A simple computation 
shows that 

n c = RlR-^ c = 2t/i£±% ( -) e~\ (2.4) 



TTj(s — 2) \S 

i.e., it agrees with what is implied by the scalings (|1.9|) and (|1.15[) . 

If f2 3> all the expressions simplify and, by using the Taylor expansion 

W(x) = |(s - 2)(1 - xf + O (|1 - x\ 3 ) , (2.5) 

and the fact that pJ F '(x; n / out ) = 0, one obtains 

PE F = H s ~ 2 )" 2 C 1 - ^n/out) 2 (1 + OQ1- x in/out |)) , (2.6) 



so that the normalization condition yields 

- 3 1 1/3 
_2tt7(s- 2)_ 

and p,Q F = C(e" 4 / 3 ri 2 / 3 ), which in turn implies that 

„TF|| 



1 ^in/out | 



(£fi)~ 2 / 3 (l + ©((eft)- 2 / 3 )) , (2.7) 



\Pn \\ L ~ im = 0((£") 2/3 ). (2.8) 
Note also that for il ^> e" 1 the density can be written as 

pZ F (x) = H s ~ 2 ) e2 ^ 2 " ^n/out) 2 - (1 - xf + O ((efi)" 1 )] + • (2.9) 
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For further convenience we introduce another reference density, which is obtained by minimizing a 
GP-like energy obtained by adding some radial kinetic energy to £ TF , i.e., if | loge| <C UJ < lu c , 

4 GP [/]:= / dx{i|V/| 2 ->W + £ - 2 [(^ + / 2 )/ 2 ]}^ (2-10) 

JR 2 

and, if fl c < Q < e~ 4 , 

£g P [f}:= [ dx{i|V/| 2 + 7 ^(x)/ 2 + £ - 2 / 4 }, (2.H) 

JR 2 



where / G ^ GP and 
We also set 



:={/GF 1 (M 2 ),/>0: ||/|| 2 = 1}. (2.12) 



E2 P := inf £S P if], E^:= inf £<f[/], (2.13) 

/eS> GP /e@ GP 



and in both cases we denote (with some abuse of notation) by go the associated unique nonnegativc 
minimizer. 

The energy E^ p (resp. E^ p ) is trivially bounded from below by E^ F (rcsp. Eq F ). Next proposition 
proves that it is also bounded from above by the same quantity up to smaller order corrections: 

Proposition 2.1 (Estimates for E^ P and E^ p ). 

If < uj < lu c as s —t then 

E™<ES P <E™+0(\\oge\), (2.14) 

while, if il c < fl <C e~ 4 , 

EE F < E% p < Eq F + 0((eil)^\ log(s 4 0)|). (2.15) 
Remark 2.1 {Radial kinetic energy of pJ F ) 

Using (|2.7p as well as (|2.8p and the Taylor expansion (|2.5|) . one easily obtains that 

El F = 0(e- 4 ^ s ), (2.16) 

which is much larger than the error term in (|2.15p . i.e., the radial kinetic energy of the (regularized) TF 
profile Pji F , for any il <C e~ 4 , since 

(e 4 n) 2 / 3 |log(e 4 r!)| « 1. 

We also remark that employing the variational equation for go one can show that Eq > Eq F + 
c(erj) 4 / 3 | log(e 4 ft)| with some c> 0. 

Proof. The result is proven by evaluating the functional £ GP (resp. £q P ) on a suitable regularizatiorj 1 " 1 ! 

of \J p^ F (resp. Pn F ) like in |CYj . The error terms in (|2.14[) and (|2.15|) are precisely given by the 
kinetic energy of such functions. For instance, if fi > fi c , it is clear that the kinetic energy of the TF 
profile far from x in / out is of order (eSl) 4 ' 3 , but one has to regularize the density in a neighborhood of 
Xi n and x ou t' By a simple optimization, one can show that the width of such regions has to be of order 
(£ 4 f2) a (ef2)~ 2 / 3 <C (a; out — £in) for some a > 0, which produces the additional factor | log(e 4 S!)| (see, e.g., 
[CPRY11 Proposition 2.1]). □ 



The regularization is necessary since the kinetic energy density is otherwise not integrable around x = x in / out . 
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With some slight abuse of notation we denote by go the minimizer of either (|2.10[) or (|2.1ip . It has 
other useful properties which can be easily proven: It is radial and strictly positive and has a unique 
maximum. Moreover one can show that it is suitably close to p TF : A trivial consequence of (|2.14p and 
H2J5H is (see, e.g., [CRYl Proposition 2.1]): 

||.9o - /°?2 F ||2 ^ C(£ 2 | loge|), if 0<c< Wc ; (2.17) 



12 TF I 

\9o ~ Pn 



2 



< O(e 10/3 n i / 3 \log(e i n)\), if f> c <f><£- 4 . (2.18) 



One can also show as in |CRY[ Proposition 2.1] that 

I, ..2 < /||p2 P L(1 + o(1)) = 0(1) j ifO<c< Wc ; 

llflblli-^) - j p TF (1)(1 + o(1)) = 0(e 2/3 fi 2/ 3) if ^ < fi <<; £ -4. ^- iJ i 

The estimates (|2.17[) and (|2.18p can actually be improved: 

Proposition 2.2 (Pointwise estimate for go). 

IfO<iu< a;J3 as £ ->• 0, #ien 

| ff0 2 W-^ F W| ^Odloger 1 ), (2.20) 

uniformly in x £ m't/i 

i:={x£l 2 : < x < x out - | logel" 1 } . (2.21) 
If on the other hand fl c < O <C £ -4 

| ffo 2 (z) - ftf < O^fiVaiJog^njia/a) ^tf^ = o(1) p TF { ^ (2 22) 

uniformly mx£ wii/i 

Aj := (xel 2 : x in + (efi)- 2 / 3 |log(e 4 r!)|- 1 < x < aw - (sQ)- 2 ^\ log^O)!" 1 } . (2.23) 

Remark 2.2 [Bulk of the condensate) 

The annulus An contains automatically the bulk of the mass provided 51 <C £~ 4 : Using (|2.7[) it is simple 
to verify that 

\\pI F \\l^a s1 ) = 1 - ^ ^^ 2 ) ( 2 - 24 ) 
pl F {x) > C{enf' 3 \ log^O)!- 1 > 0, (2.25) 



and 



for any x G An- The same obviously holds true for A u as well: 

^ F (x)>C| logel" 1 >0, Hp^Li^) =l-0(|log £ |- 2 ). (2.26) 

As a direct consequence of the pointwise estimates (|2.20[) and (|2.22[) , A u and An contain the bulk of the 
mass of go too: 

IkllL^) = 1-O(|loge|- 1 ), if 0<uj<u c , (2.27) 
!l3o|ll 2( ^) = l-0(|log( e 4 fi)r 2 ), if n c <n«s-\ (2.28) 



12 The condition o> < u) c should be thought of for u) = 0(e 1 ) as the stronger requirement u> = ljqe 1 for some constant 
ljo < u) c £. Similarly the condition f2 > Q c means actually O > (1 — o(l))O c and, if f2 is below f! c , x ln has to be set equal 
to 0. 
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Proof. The proof is done exactly as in |CRYl Proposition 2.6] (see also jCPRYll Proposition 2.3]), i.e., by 
identifying local super- and subsolutions to the variational equation solved by go, which can be rewritten 

- g'l - x- l g' = As' 2 [p - g 2 } g , (2.29) 

where 

5(x) ._ / \ W -* S + , if < u < Uc , 

' 1 3 [e 2 £g P ~ je 2 n 2 W(x)] , if O c < n « e-\ 1 ' 

and the chemical potentials and /1q P are fixed by the L 2 -normalization of go. Note that p differs 
from the TF densities pJ F and p]^ F only by the chemical potential. However combining (|2.17[) . (|2.18l) 
with (|2.19|) . one can easily obtain the following estimate of the difference between the chemical potentials: 

|AS P -^ F | ^(e^llogel 1 / 2 ), if 0<u<u c ; (2.31) 
|A8 P -Mn F | < O(^|log(e 4 0)| 1/2 ), if n c <ft<£- 4 . (2.32) 

This in turn yields a pointwise estimate of the difference between p and the TF densities, which guarantees 
that p is strictly positive inside A u or An- 

The rest of the proof consists in a local analysis of the variational equation (|2.29p on the scales e or 
e 2 / 3 $l _1 / 3 | log(e 4 f2)|, if lo < uj c or Q c < <C e -4 respectively. Note that in this second case the blow-up 
scale remains much smaller than (ef2)~ 2 / 3 | log(e 4 f2)| _1 which appears in the definition of An- □ 

We conclude this section by stating another pointwise estimate of go which shows the exponential 
smallncss in e of go outside of supp (p TF ) • 

Proposition 2.3 (Exponential smallness of go). 

There exists a constant C > such that, if | log e| <C u> and 11 <C e -4 as e — > 0, 

2 fCcxp{-£- 1 / 2 dist(x,supp(^ F ))}, if [ log e] < tu < UJ C , 

9o{X> ~ }C(rf}) 2 / 3 exp {-J2V2 dist(x,supp(pTF))} ; if 0c < n « e" 4 , 1 ' 

for any x G R 2 . 

Proof. See, e.g., |CRY[ Propostion 2.2]: For instance in the first case, i.e., when | loge| <C w < w C5 on e can 
show that if x > x ou t + y/e, the r.h.s. of (|2.29|) is negative and C exp{— e _1 / 2 (x — x ou t)} is a supersolution 
to (|2.29p for x > x out and a suitably large constant C > 0. In the other case when il c < <C £ -4 one can 
investigate (|2.29p separately for x < a;; n and a; > x ou t and identify appropriate supcrsolutions there. □ 



2.2 Energy Asymptotics and Emergence of a Hole 

In this section we first prove Thcorcm ll.il i.e., the asymptotics of the GP ground state energy, and show 
how this implies a phase transition at f2 C2 , that is the occurrence of a macroscopic hole in the condensate. 

Proof of Theorem \l.ll For the sake of brevity we omit the proof of (| 1 . 3T[) and focus on the regime 
Sl c < f2 <§C e~ 4 . The strategy as well as most of the estimates are basically the same for (|1.37[) and (|1.38D 
but the vanishing of the TF density for x < Xi n makes the second case slightly tougher to deal with. 

As usual the result follows from a comparison of appropriate upper and lower bounds to the GP 
energy E^ p , but in order to simplify the analysis it is convenient to extract the leading order term 
from the outset. This can be done by using a very standard trick (sec, e.g., LM ), i.e., a splitting of the 
energy: Exploiting the positivity of go, we define a function uq G L 2 oc (R 2 ) as 



i>n ( x ) = : 9o{x)u (x), 



(2.34) 
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and compute, using the variational equation (|2.29p solved by go, 

E$ F = E§ f +F[uq], (2.35) 
where T is the weighted GL-type functional 

/ dx 5o 2 {i |(V - iAn) u| 2 + e- 2 5o 2 (1 - M 2 ) 2 } . (2.36) 



T\u\ := 



We can then focus on the reduced energy .F[«o] and prove upper and lower bounds for it. Note that in 
the upper bound there is a normalization condition to fulfill, i.e., 

dx ff g|uo| 2 = l. (2.37) 

In the lower bound on the opposite we are going to remove this restriction and minimize w.r.t. any 
u g L 2 (An), where An is given by (|2"25|) . 

Upper bound: The trial function which is going to be used in the upper bound proof is very similar 
to the one in |CY1 Proof of Proposition 4.1]: Roughly speaking we pick a function Utrial which is approx- 
imately a phase factor, i.e., such that |utriai| — 1 5 but at the same time contains vortices of unit degree 
on a regular lattice intersected with the support of the TF density p FF . More precisely we denote by £ a 
regular lattice with fundamental cell Q (triangular, square or hexagonal) and lattice spacing £ and define 

"triai(x) := c £(x) u(x), (2.38) 

where c ~ 1 is a normalization constant to enforce (|2.37[) . v a phase factor (here we denote a point 
x = (xi, X2) £ R 2 by the complex notation Q = x\ + ix%) given by 

V W ;= II =:exp{^(x)}, (2.39) 

and ^ regularizes the singularities of v at the lattice points 



,1, if |x - Xj| > t, for all x, g C, 

U-^x-xd, if Ix-Xil <t. 



Here t > is a variational parameter (vortex core radius), such that 

t<fT 1/2 . (2.41) 
We also choose the lattice spacing £ in such a way that 

|Q|=7rn-\ (2.42) 

which guarantees the 'neutrality' of each cell of the lattice in the 'electrostatic analogy' discussed in |CY1 
Proof of Proposition 4.1]. Since |supp (pJi F ) \ = 0((efl)~ 2 ^ 3 ), the total number of lattice points inside 
supp (/Of2 F ) is 0(e~ 2 / 3 S! -1 / 6 ) ^> 1, as long as <C e~ 4 , and it is not difficult to verify that 

c = 1 + 0(e- 2/3 n-^ 6 t 2 ) = 1 + o(l), (2.43) 

thanks to (|2Hj) . 

The evaluation of ^[utrial] can be decomposed into two different estimates of the kinetic energy and 
interaction energy respectively. We are going to show that the former satisfies the inequality 

|/ dx.g 2 |(V-zA^) Utrial | 2 < i(l + o(l))fi|log(* 2 n)|. (2.44) 

Jm 2 



Critical Speeds in the GP Theory - CPRY - March 27, 2012 



21 



A simple computation yields (recall that £ is real) 

§/ dx 5 2|( V -zAoK r iai| 2 = ic 2 f dx ff0 2 £ 2 |V^-A o | 2 + c 2 V f dx<7 2 |Ve| 2 < 
Jr 2 Jr 2 x7ec^l x ^ x 'l^* 

±(l + o(l)) / dx 5 2 £ 2 |V<^- A n \ 2 + Cn^ 2 t + 0(e°°), (2.45) 



where .A is chosen in such a way that (|2.33[) yields the pointwisc bound 43 ! go — 0(s°°) outside A. A 
possible choice is 

i:={x€K 2 : x in - fr 1 / 2 ! log(e 4 fi)| 2 < x < x out + fi" 1 / 2 ! log(e 4 fi)| 2 } . (2.46) 

Now we can apply the electrostatic computation in |CYj to estimate the leading term in (|2.45[) : The first 
term on the r.h.s. of (|2.45l) becomes 

| / dx 5o 2 £ 2 |V^-A u | 2 <±(l + 0(ifi 1/2 )) V sup g 2 (x) \ir\A\\\og(t 2 n)\+0(l)} , (2.47) 

as in |CY1 Eq. (4.37)]. In the above expression we have actually overestimated the contribution in 
A \ An since v contains no vortices there and therefore one would easily get an upper bound without the 
logarithmic term in the cell contained in A \ An- 

Using a Riemann sum approximation for the integral of p 2 in A we can replace the sum in (|2.47|) and 
obtain 

i / dx 5 2 £ 2 |V^A | 2 < i(l + O (l))ft[|log(i 2 ft)|+0(l)] , (2.48) 

thanks to (|2.4ip and 

ho\\l H A) ^ WpoWlHAo) + °M + yo\\l H A\A n ) < 1 + °(!) + 0(| log(£ 4 ri)|- 2 ) = 1 + o(l), (2.49) 

since g 2 < (eft) 2 / 3 | log(e 4 fi)|- 1 in A\Aq and n-V2| log( £ 4 Q)| 2 < (e«)- 2 / 3 | log(e 4 fi)|- 1 for any (1 < e~ 4 . 

The upper bound (|2.44[) is then proven and it remains to compute the interaction energy of ittriai : For 
any fi c < O -c £~ 4 one can show that 



e- 2 



j2 dx gl (l - K, al | 2 ) < Ce-^n b 'h 2 . (2.50) 
Indeed we note that |ittriai| = 1 in R 2 \ An, so that 



/ dx.g 4 (l - Kiail 2 ) = J2 [ 

Jr 2 v y „JtZa-Jb 



,2^ 



dx ff S l-^lx-Xil^ < 

— An JB( Xi ,t) V 7 

C( e fi) 4 / 3 t 2 |Q|-Vn| < Ce 2 ^ 3 t 2 . (2.51) 

Combining (j2~35|) with ([2lH|) and (|2~5U)l . we finally obtain 

Eg P <Eg p + ±(l + o(l)) | log (t 2 ft) | + Ce-^rf'h 2 + o(l) < 

+ 5 (! + °(!)) "I lo g(^)l + <W ( 2 - 52 ) 



13 Note that 112. 33 II also implies an exponential decay in x of go, which is used in order to estimate the remainder in 12.451 1 
14 The symbol C(e°°) denotes a quantity vanishing faster than any power of e. 
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where we have chosen the vortex core radius t = £ 2 / 3 fi -1 / 3 and used (|2.15p . Note that our choice of t is 
compatible with (|2.41[) as long as O -c £~ 4 - Moreover the radial kinetic energy of the profile go, which is 
0((efl) 4 ^ 3 \ log(£ 4 n)|) according to (|2.15p . is much smaller than the remainder 0(fl) if again fl <C £ -4 . 

Lower Bound: A lower estimate matching (|2.52[) can be proven by following the same strategy as in 
jCYl Section 5], i.e., first decomposing the support of into cells, then blowing-up the energy .F[uo] on 
each cell and finally using a lower bound for the GL functional with an applied magnetic field between 
the second and third critical values (see, e.g., |SS1[ ISS2] ). 

The starting point is a restriction of the integration in J- to the domain An, which exploits the 
positivity of the integrand. Then we decompose An into cells: Let £ be the lattice 15 ! 



C = 



jxj = (ml, n£), m, n e Z j (2.53) 



with cells Qi and lattice spacing satisfying 

Q^l log(e 4 f!)| < P < £ - 4/3 fi- 4/3 , (2.54) 

we easily obtain the lower bound 



F[u ]> £ J- dx^{i|(V-iA n ) U o| 2 +£- 2 3o 2 (l-K| 2 ) 2 }> 



Xi.Qi CAn 



(i-o(i)) J2 sSW^N. ( 2 - 55 ) 

where we have used (|2.22p and (|2.25|) . which imply the lower bound 

gl(x) > C(efi) 2/3 | log^O)!" 1 , for any x € An, (2.56) 

and introduced the functional 

•F (i) M:= J dx{i|(V- i Ao) U o| 2 + e- 4/3 ^ 2/3 |log(£ 4 «)r 1 (l-|«o| 2 ) 2 }. (2.57) 

Note that the conditions (|2.54[) on the spacing £, which are essential for the proof, are compatible as long 
as il <C £~ 4 : The upper bound on £ is needed in order that the spacing is much smaller than the width of 
supp (pj^) as required for the Riemann approximation we shall use. On the other hand, the lower bound 
in (|2.54p implies that the cells of the lattice C are much larger than those of the lattice considered in the 
proof of the upper bound. The reason is that we would like each cell of the new lattice C to contain a 
large number of vortices, diverging to oo as e — > 0. In fact, this is crucial in order to make the whole 
proof strategy work since the GL estimate from jSSll ISS2] that we are going to use holds only under this 
condition. 

We now blow-up the energy in the cell Qi and set s := £^ 1 (x — x,), 

t2 (s) := w ( Xl + Ix), A(s) := £ A fi (x 4 + &), (2.58) 

so that becomes 

T^[u Q ]=T^[u Q ]:= ds ji|(v-zA)u | 2 + e - 2 (l-K| 2 ) 2 }, (2-59) 



°For convenience we have chosen a square regular lattice but any regular lattice covering the space would have done the 
job. 
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where Q\ stands for a unit square centered at the origin and we have introduced a new small parameter 
e defined as 

e := th^n-^l \og(e 4 n)\^ 2 , (2.60) 

which is much less than (e 4 ^) 1 / 3 ! \og(e 4 fl)\ 1 ^ 2 <C 1 if SI <C £~ 4 thanks to (|2.54l) . The rescaled vector 
potential A is explicitly given by 

A(s)=We,Ax 1 +ft?e 2 As, (2.61) 
and the corresponding magnetic field is 

h cx = curl A = Q£ 2 . (2.62) 
Employing gauge invariance one obtains the lower bound 



P l) [u ]> inf / ds|i|(V-^ cx e,AsH 2 + e- 2 (f-M 2 ) 2 } 
uem(Q 1 )J Ql I 2 J 



(2.63) 



where the applied magnetic field h ex satisfies the conditions 

|loge|<h cx <e- 2 , (2.64) 

thanks to (f2~54"|) . 

We can now borrow a lower bound to the GL functional with an applied magnetic field between the 
first and second critical fields (see, e.g., [SSI , SS2 ) to estimate the r.h.s. of (|2.63D and obtain 

•F W N > (1 -o(l))h cx log (e" 1 /*- 1 / 2 ) > i(l- (l))fW 2 |log( £ 4 ft)|. (2.65) 

Putting together the above estimate with (|2.55[) and using again the Riemann sum approximation we 
finally get 

•FN > i(l - o(l))fi| log( £ 4 Sl)| ||5o|| 2 L2( ^) > i(l - o(l))fi| log( £ 4 fi)|, (2.66) 

by (|2.28p . The splitting (|2.35p as well as the trivial bound -E GP > E^ F hence complete the proof of the 
lower bound. □ 

The occurrence of a hole in the bulk of the condensate for f2 > Sl c now follows as a straightforward 
consequence of the energy asymptotics proven in Theorem ll.il 



Proof of Proposition 11.11 The proof is very similar to the proof of Proposition 12.31 (see also jCRYl Pro- 
postion 2.2]): By setting p GP := |"0n P | 2 one easily obtains from (|1.25[) the inequality 

- ±Ap GP < 2e- 2 [p GP - p GP ] p GP , (2.67) 

where 

p GP (x):=I[eVH P -jeWW(x)]. (2.68) 
Now the energy asymptotics (|1.38|) implies exactly as in the derivation of (|2.18|) and (|2.19D that 

|^ p | < C(sil) 2 /\ ||p GP - pji F \\l < Cs 2 n\ log(s 4 0)|, (2.69) 

which in turn (see, e.g., jCRYl Proposition 2.1]) yield the estimate 

ImH P -^ F | < Ce-^^llog^ 4 ^)! 1 / 2 . (2.70) 
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Now if > f2 c , one can pick some < S <C 1 so that the above inequality implies (recall (12.51) ) 

p GP (x) < Ce 4 ^ e \ log(e 4 ft)r/ 2 - C{en) 2 'H\ (2.71) 

for any < x < x- m — 8(eQ)~ 2 ^ 3 . Hence if one takes for instance 

S ^ (e 4 n)^ 12 \\og(e 4 n)\ (2.72) 

that is <C 1 as long as f2 <C £ -4 , the inequality (|2.67|) becomes 

- ^Ap GP + Cs- 2/3 n 5/e p GP < 0, for any < x < x in ~ 6(en)- 2/3 . (2.73) 

To complete the proof it suffices therefore to note that the function Ccxp{— e^ 1 ^ 3 fl 5 ^ 12 (xi n — x)} is a 
supersolution to the equation —Af + Ce -2 / 3 5 / 6 / = in the ball < x < X[ n with boundary condition 
/ = C on dB Xin . One thus obtains 

p GP (x)<C\\p GP \\ oo e^{-e- 1 ^/ 12 (x in -x)}, (2.74) 

for < x < Xi B but, if Q > HqE^ 1 for Qq > efl c , then x ln > C > and one can always find some 
< Xhoic < ^in such that the above upper bound yields the desired estimate. □ 

2.3 Uniform Distribution of Vorticity 



The energy asymptotics proven in Theorem 11.11 is the main ingredient in the proof of the uniform distri- 
bution of vorticity: 

Proof of Theorem \1.2i The proof follows very closely jCPRYli Proof of Theorem 1.1], which in turn relies 
on [SS11 Proposition 5.1]. For the sake of brevity we only enlighten the major differences and assume 
that fl c < < e~ 4 . 

The starting point is the combination of the cell decomposition (|2.55[) with the global upper bound 
()2.52p . which gives (recall that £ is a regular square lattice with spacing I satisfying (|2.54[1 ) 

9 2 o(x t )^[u ]<^l + v ) 5o(^K 2 ^|log(e 4 ^)l, (2.75) 

for some quantity r\ := 7y(e, fi) such that rj — > as e — > 0, f2 <C £~ 4 . 

In order to localize the above bound, one needs some control on the density go and to this purpose 
we set 

(3:=f3(e,n)=m a x(\\ogr ] \-\ llog^O)!" 1 ), (2.76) 

and define 7\Lbuik as in (jl.42p . i.e., 7?.buik = {x : x- m + (efl)~ 2 / 3 (3 < x < x out — (£^)~ 2 ^ 3 /?}, so that 
Hbuik C An and 

g 2 (x) > C(eft) 2/3 /3, for any x e ^ bulk . (2.77) 

Such a bound unfortunately is not sufficient to complete the localization of the energy and a certain 
amount of bad cells has to be rejected: We say that a cell Qi is good if 

F (l) [u ]<\ 1 + min (^/ij, e 2 VU^j nl 2 \ log(e 4 fi)|. (2.78) 

Now given any set S C 7?-buik with area much larger than I 2 (|2.75[) implies that the total number of bad 
cells contained inside S (denoted by Nb) is much smaller than the total number of cells N, i.e., more 
precisely 

N B < Vvfi^N < min L/^log??!, e 2 VTi\ log(e 4 rj)|] N < N. (2.79) 

On the other hand one can use the upper bound ()2.78[) inside good cells to prove the final result exactly 
as in jCPRYli Proof of Theorem 1.1]. Note that the choice I = Q- 1 ^ 2 \ log(e 4 ft)| is compatible with ([234"]) 
and yields the lower bound on the measure of S, i.e., \S\ ^> £ 2 = f2 _1 | log(e 4 f£)| 2 . □ 
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3 The Giant Vortex Regime Q ~ e 



-4 



Since throughout this section we deal only with the regime fi ~ e 4 , we remove any suffix from tf) GF , 
£^f p , etc. for the sake of simplicity and denote those quantities by ip GP , E GP etc.. 

3.1 Preliminary Estimates for E GP and ip GP 

We first introduce and recall some notation: As in (|1.53[) we define a parameter a as 

a 2 :=4 + 7(s-2), 

and a vector field B by 

B = (fir - |fijr _1 ) e#. (3.1) 
A Taylor expansion of the potential U(x) (see (|1.50[1 ) around x = 1 (e.g., for i < x < |) yields 



U(x) = \B 2 (1) + (1 - LfiJ 2 fi" 2 ) (x - 1) + \ [1 + 3LfiJ 2 fT 2 + 7(* " 2)] (x - l) 2 + O ((a; - l) 3 ) 

i 

i 



1 a 2 (x - l) 2 + 0(n- 2 ) + O ( {x - if) . (3.2) 



Moreover we denote by g osc the normalized ground state of the one-dimensional harmonic oscillator 

(yen) 

#osc := -\d 2 y + \a 2 y 2 , (3.3) 

5osc(y) ^vr-VV^expj-iay 2 }. (3.4) 
The starting point is a rough but useful upper bound on E which has several important consequences: 

Proposition 3.1 (Estimates for E GP and ijj GP ). 
If CI = fio£ -4 with fio > and e is sufficiently small, 

E GP <0(fl), W^Wl <0(e 2 n). (3.5) 

Moreover there exists a finite constant C such that, for any x £ M. 2 , 

|V GP (x)| 2 < Ce^expj-fi^ll-a:!}. (3.6) 

Proof. In order to prove the upper bound on E GP , it suffices to test the GP functional on a suitable 
rescaling of the harmonic oscillator ground state (|3.4p times a phase containing a giant vortex of degree 
|_fij at the origin. More precisely we set 

^triai(x) :=./triai(z)e W , (3.7) 

fn 1 /* 5oBC (nVa(i_ a .)) ) if a; > 1 - e, 
/triai(^) := c n I f{x), if 1 - 2e < x < 1 - s, (3.8) 

[o, if0<x<l-2e, 

where is a normalization constant and / a smooth function such that 

/(l - 2e) = 0, /(l -e) = ffosc(fi 1/2 £) = 0(e°°), 
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which also allows to assume that ||V/|| = 0(e°°). A simple computation gives 

/>oo 

1 = HVtriaiHa = 2^4^/2 / dxxg^n^ix- l)) + 0(e°°) = 

Jl-E 

dy(l + Q.- l / 2 y)gi ac {y) + 0{e°°) = 2*4 /°° dy g 2 osc (y) + ©(IT 1 / 2 ), 



2-KCt 



so that 



(27T)- 1 (l+O^" 1 / 2 )) 



(3.9) 



By using p. 21) and the exponential decay of (fosc one has 

poo 

£8 P Atrial] = 2*4 / dx .T {i (/ t ' rial ) 2 + « 2 C/(x)/ t 2 rial + 4^ 2 /t 4 rial} 



2tt4 



dx Z {l (/ t ' rial ) 2 + fi 2 f/(x)/ 2 rial + 4£" 2 /t 4 rial} + 



2*4 n 



-efiV 



1/2 I. 



+ 4^ 



osc T" H2 



2*4° 



a 1 



-1/2n 



(3.10) 



2 2?r V 2*^0 

The second statement in (|3.5D is a consequence of the variational equation (|1.25[) which yields 

- ±Ap GP < e- 2 [e 2 M GP - je 2 Q 2 W - 2p GP ] p GP < e~ 2 [s 2 p GP - 2p GP ] p GP , (3.11) 

where 

(3.12) 



P 



GP L/.GP 



and we have used the positivity of W. Now since Ap GP < at any maximum point of p , it immediately 
follows that 

I|P GP L * ^ GP , (3.13) 
which together with the definition of the chemical potential implies 



M GP = E GP + e -a ||^GPir < E GP + 1 M GP ; 



so that by ([3TTU|) 

/i GP < 2£ GP < (3.14) 

and thus Q3.5P is proven. 

In order to prove the pointwisc estimate (|3.6[) we exhibit an explicit supersolution to the variational 
equation (|1.25[) . We first consider a ball B(l — as 2 ) centered at the origin of radius 1 — as 2 . with a = 0(1) 
sufficiently large: The monotonicity of W together with the Taylor expansion W(x) = |(s — 2)(1 — x) 2 + 
O ((1 - a;) 3 ), yields the bound 

W(x) > Ca 2 e 4 , 

for some C > and for any < x < 1 — as 2 . Therefore inside B(l — as 2 ) one has by (|3.11[) and (|3.14[) 

-iA P GP < e~ 2 [eV GP - Ca 2 e 6 n 2 } p GP < -\Ca 2 e^ 2 P GP , 
if a is taken sufficiently large. In the same region the function 



/» P (x) := C a Wp^W^ exp {-^(l - x)) , 



(3.15) 
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satisfies 



h^fsup + hCa 2 £^fl 2 /sup — h 



-n 1/2 -n x + c*aVr2 : 



/sup > o, 



if a is again large enough. Therefore / sup provides a supersolution, since at the boundary 923(1 — ae 2 ) 
one has 



/ sup (l - as 2 ) = C a Wp^W^ exp {-a^ /2 } 



= P 



GPI 



if C a is taken equal to exp{afty 2 }. In conclusion, by the maximum principle (see, e.g., [TJ Exercise 2, 
p. 317] or P Theorem 1, p. 508]) p GP (a;) < / sup (:r) for any x e B^ ae 2 but the result can be trivially 
extended to the whole ball B\ thanks to the monotonicity of / SU p- 

To complete the proof of the pointwise estimate (|3.6p . one has simply to repeat the argument above 
in the complement of the region Z3 1+ae 2, exploiting the monotonicity of W there, and replace the super- 
solution / SU p with the function 



a||p GP Lexp{-^ (a; -i)} 



□ 



The decay of ip GP outside of the bulk of the condensate can be estimated more precisely by exploiting 
the Taylor expansion (|2.5|) : 

Proposition 3.2 (Exponential decay of ip GP ). 

If O = £IqE ~ 4 with Ho > and 1 <C A <C | loge| b for some b > as e — > 0, there exist two finite constants 
c > and C < oo independent of e and X such that 



|V> GP W| 2 < Ce- 2 exp {-eft (|1 -x\- e 2 A) 2 } , 
for any x satisfying the condition 

e 2 X < |1 - x\ < e 2 X a , with a > 2. 

Proof. We first consider the region x < 1 — e 2 A and prove that the function 



fsup(x) 



„GP| 



exp {-cO(l - e 2 A - a;) 2 } + exp {-^/n^ l X 2a+1 } 



(3.16) 



(3.17) 



(3.18) 



provides a supersolution to (|3.11[) in the region e 2 A < 1 — x < e 2 X 2a+1 , for any a > 2: 

[-|A + 7 ft 2 M/(x) + 2 £ - 2 / 2 up ( a; ) - p GP ] / sup > 

[|cfi + |(s - 2 - c 2 )ft 2 (l - x) 2 - Ce- 2 A 6Q+3 - CO] /s U p > 



--4 r i 



(s - 2 - c 2 )ft 2 A 2 - Cfto] /su P > 0, (3.19) 



if c 2 < s — 2 and e is sufficiently small (recall that by assumption A — > oo as e — > 0). Moreover at the 
boundary 



/su P (l - £ 2 A) > lip 



GPI 



> P ^(l-e 2 X), 



f^p{l~s 2 \ 2a+l )>p w (l-E 2 X 2 ^) 



(3.20) 



thanks to 



Therefore p GP < / sup and in the region (|3.17p 

exp |cO(l - e 2 X - x) 2 - V^A 2 ^ 1 } < exp ^-^/Ur]x 2a {X - C)\ = o(l). 



The proof for x > 1 + e 2 X is identical. 



□ 
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3.2 The Giant Vortex Density Profile 



In this section we investigate the properties of the giant vortex profile and the associated energy functional 
defined in ([05)1 . 

Proposition 3.3 (Minimization of £ gv ). 

There exists a minimizer g gv of p. 491) that is unique up to a sign, radial and can be chosen to be strictly 
positive for x^O. It solves the variational equation 



\Ag gv + ti 1 U(x)g gv + 2e 2 5 f v = /i gv 5gv , 



(3.21) 



with /i gv = E sv + e 2 1 1 £?gv 1 1 4 - Moreover g gv has a unique maximum and, if VI = floe 4 with Qq > 0, 



In addition there exists a finite constant C such that 



sup 

iei 2 



27T11 



- 1/4 Sgv(*)-3o S c (^(l-a-)) 



< CCl, 



-1/4-1 



(3.22) 



(3.23) 



Finally g gv decays exponentially far from x = 1: // 1 < A < | loge| h for some b > there exist constants 
c > and C < oo independent of e and A such that, for any e 2 X < |1 — x\ < e 2 X a , a>2, 



5gvO*0 < C*e- 2 exp{-cfi(|l-2;| -e 2 A) 2 } 



(3.24) 



Proof. All the properties of g gv can be deduced by standard arguments (see, e.g., |CRY1 Proposition 2.3] 
or [CPRYll Proposition 4.1]). For instance the existence of a unique maximum point of g gv can be proven 
by an adaptation of the argument in |CPRYli Proposition 2.2], while the pointwise estimate (|3.24[) can 
be proven exactly as (|3.16l) . 

The energy upper bound in (|3.22p . i.e., 



i 



2 2?r V 27rfto 



-1/2n 



(3.25) 



has already been proven in (|3.10[) . In order to prove a lower bound matching with (|3.10p , one needs first 
to show that 



9gv 



(x) (2n)-^n^g osc - *)) J] = O(fi - 1/4 ) 



(3.26) 



After a rescaling of all lengths in the functional (|1.49[) . i.e., setting y := S7 1 / 2 (1 — x), and denoting 



g gv (y) := ^tt) 1 / 2 ^ ^ (l - fl^y) , 



(3.27) 



one obtains 



l l/2 



d V {l l V 5gv| 2 + \a 2 y 2 g 2 gv 



+ (27T)"^- 1/ \ 4 V }> 
Q, (l - O^- 1 ' 2 )) (^1 H osc \g gv ) (3.28) 



where we have used the Taylor expansion (|3.2[) and H osc is defined in (|3.3|) . Strictly speaking the 
expectation value on the r.h.s. of the above expressions is computed in L 2 {— oo, ^fl 1 ^ 2 ) but, exploiting 
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the exponential smallness (|3.24[) of g gv at y = i£l 1//2 , it is not difficult to slightly extend g gv smoothly to 
y larger than \^t 1 ^ 2 in such a way that the mean value can be computed in L 2 (R). The error 0(e°°) due 

— 1/2 

to this procedure can safely be included in the remainder O($l ). 

We now decompose g gv in a Fourier series in terms of the normalized cigenfunctions g n {y), n E N, of 
the harmonic oscillator (|3.3[) (with go = g osc ), i.e., set 

oo 

9gv(v) = a n9n(y), (3.29) 

n=0 

and use the upper bounds above to get 

^ + yEW 2 ^E("+^) = «W Iflkr) < f + f 7^ + 0(^ 1/2 ), (3.30) 

n=l n=0 71 ' 71-0 

which immediately implies 

oo 

^K| 2 <CQ„ 1/2 , (3.31) 

n=l 

since £ |a„| 2 = 1 - 0(e°°) by ([3T21]) . Therefore (|3~2"6"j) is proven and, since HffgvIL < CfiJ /2 (see below), 
one also has 

Hffgvlll > llflosc||4 - Cnl /2 \\g gv - g osc \\, 2 - 0(e°°) > \\g osc \\ 4 - C^ 4 , (3.32) 
which yields (j3~2"2"j) . 

In fact (|3.3ip has much stronger consequences than (|3.26p and, if we define for any g £ L 2 (R) 

NlL : =<<?l#o S c|<?>, (3.33) 

llffgv - 5osc|| 2 sc = (ffgvl H osc \g gv ) - a (g gv \g osc ) + \a < Cn~ 1/2 , (3.34) 
where we have used (|3.22[) and the inequality 

(5gv|5osc) > 1 - Cfi(7 1/2 , 

which is a trivial consequence of Q3.26p . On the other hand the Sobolev embedding theorem in R 
immediately implies 

llfflli-(K) <C\\g\\ H i {R) < C||.g|| osc , 

and therefore Q3.23P follows from (|3.34p . Note that the prefactor in (|3.23p is bounded independently of 
Clo because it behaves as fi 1 ^ 4 f2 1 / 4 . □ 

For technical reasons which will be clearer later we also consider a functional with a different integra- 
tion domain, i.e., 

A Ve :={x£l 2 : 1 - e% <x<l + e 2 r] £ } , (3.35) 

where 

|log £ |« % «|loge| b (3.36) 

with b > 1 is a parameter which will be fixed later as rj E = | loge| 3 / 2 . Note that since ?y £ ^> 1 the domain 
A„ e expands on a scale e~ 2 , which ensures that it contains the bulk of the mass. We define 

9% := {/ 6 H 1 (A, h ) :/ = /*, ||/|| i2(A , e) = l} (3.37) 



one has 
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and set, for any / £ 

:= / dx{i |V./f + tfU(x)f +£- 2 / 4 } . (3.38) 

Wc recall that 

U(x) = \B 2 {x)+ 1 W{x) = i|fij 2 ^" 2 ar 2 + ix 2 + i 7 (x s - 1) - ±7 (x 2 - l) - > 0. 

The ground state energy is 

Ef := inf £f [/] (3.39) 

and we denote by g, h any associated minimizer. 
Proposition 3.4 (Minimization of £|J). 

There exists a minimizer g Vs of (|3.39[) i/ia£ is unique up to a sign, radial and can be chosen to be strictly 
positive for x ^ 0. Inside A Ve it solves the variational equation 

- \ Ag, h + n 2 U(x)g„ E + 2e~ 2 gl = ^g Ve , (3.40) 

with boundary conditions <7^(l±e 2 ?7 e ) = and nf^ — E^ +e~ 2 \\g VE \\ 4 - Moreover g„ E has a unique global 
maximum at x max . If ft — f2o£ -4 with VLq > 0, then Hff^H^ = 0(e 2 £l) and 

E% = (l + 0(e°°))£; gv . (3.41) 

Proof. The only result which deserves a discussion is the energy upper bound (|3.41[) since anything else 
can be dealt with as in the proof of Propositions 13.11 and 13.31 However it suffices to test the functional 
on a trial function of the form c £ g gv , where c £ ensures the normalization in L 2 (_4,, e ), to get 

[Ce3gv] < C\E^ [ ffgv ] = c\E*\ 

since c e > 1 and the energy is positive so that we can extend the integration domain from A„ E to the 
whole of R 2 . It remains then to estimate the normalization constant, but the exponential decay (|3.24D 
together with the conditions (|3.36[) on r] e guarantees that c e = 1 + 0(e°°). □ 

In the next proposition we prove the analogue of (|3.24[) . i.e., an estimate of the decaying rate of g„ E . 
Actually we also state more refined pointwisc estimates of g VE which will be crucial in the proof of the 
absence of vortices. 

Proposition 3.5 (Pointwise estimates for g„ E ). 

IfCl = S1q£ — 4 with flo > 0, there exists constants C, c > 0, such that for any x with |1 — a:| > e 2 r)l^ 4 , 




(3.42) 



Moreover if in addition J7o > with Qq = 0(1) large enough but independent of e, then 

= (l + 0(n- 1/4 )) (2^)- 1 / 2 r» 1 / 4 , 9osc (V/ 2 (l - a;)) (3.43) 

uniformly in x £ Ayjjj. 
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Remark 3.1 {Maximum of g, h ) 

A straightforward consequence of p.43[) is the estimate 

< ^- 3/2 V^(l + c«o 1/4 ) . (3.44) 

On the other hand since g' (x max ) = and g'S(x max ) < 0, the variational equation (|3.40j) implies the 
following estimate of the maximum position 

(1 - .T max ) 2 < 2 (any 2 - 2e- 2 gl(x m a X )] < (aiiy 1 (l + Cn^ 2 ) . 

In fact, (|3.48[) below yields a better estimate: 

ln-*/ 2 Van (l - Cfio 1 / 4 ) = (l-Cn- 1 ^)) (2k)-iW g 2 osc (0) < £(1) < 

ff 2 > max ) < ivr" 3 / 2 ^ (l + CO" 174 ) exp {-aQ (1 - x mgx ) 2 } , (3.45) 

which implies 

|l-aw| < C% 5/8 e 2 . (3.46) 



Proof of Provosition Iff. 51 The proof of (|3.42|) is basically identical to the proof of (|3.24[) . The only 
difference is due to the fact that g ne does not vanish at the boundary of A Ve ■ This implies that any 
legitimate supersolution to the variational equation (|3.40[) must be larger than g Ve at dA Vc ■ Alternatively 
one can simply repeat the proof of ()3 . 1 6[) or (|3.24|> , but extend g Ve outside A Ve to a function which decays 
to zero. A simple way to achieve that is to construct a function g ric which equals g Ve in A ne and satisfies 
the inequality 

- iA^ + n 2 U(x)~ griE + 2e- 2 ~gl < (i%g Vs , (3.47) 

in the region e 2 r\ e < |1 — x\ < e 2 ri 2 , together with Dirichlet boundary conditions g Ve (1 ± e 2 r) 2 ) — 0. Note 
that g Ve is a priori only a weak subsolution because it is in general only continuous at 1 ± e 2 r\ e . However 
any supersolution to p.47[) in the interval |1 — x\ < e 2 rj 2 provides a pointwise upper bound to g T?£ in A„ e 
and, acting as in the proof of (|3.24[) . it is not difficult to realize that the r.h.s. of (|3 .42[) is a supersolution. 

The second statement is obtained by combining the application of the maximum principle to the 
variational equation (|3.40[) far from x max with elliptic estimates like (|3.23p in the central region around 
Xmax where g„ e is large enough. A direct consequence of the decay estimate (|3.42[) is indeed that one 
can apply the analysis contained in the proof of Proposition 13.31 In particular the analogues of (|3.28|) - 
(13.341) hold true, which yield as in the proof of (|3.23[) . 

I^fi-^VO*) - <w (r> 1/2 (l - a?)) 

Note that (|3.42[) is needed in order to extend the analysis to the whole real line exactly as in the discussion 
after (|3.28|) . The estimate above is already sufficient to obtain (|3.43|) in the region where g Vs is large 
enough, i.e., sufficiently close to x max . However far from x max a different approach is needed. We start 
by proving the following upper bound: For any x € A Ve such that 

\l- x\>{att r) )- 1 e 2 1 (3.49) 

there exists a constant c so that 

9 2 E (x) < (l + cO,7 1/4 ) i7r- 3 / 2 v /^exp{-afi(|l-2:|-2(afio)~ 1 £ 2 ) 2 }. (3.50) 



< Cfi, 



-1/4, 



(3.48) 
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As above we use an extension g Ve of g Ve to the interval |1 — x\ < e 2 r] 2 and restrict for simplicity the 
analysis to the region 

(3.51) 



<x< 1- (afio) 1 £ 2 , 



where we consider the function 



/ 2 up (x) := (l + dV /4 ) i 7 r- 3 / 2 y^cxp|-ar!(l-(«r!o)" 1 e 2 -.T) 2 |, (3.52) 

At the boundary of the interval (|3.51|) one certainly has g, h < f sup thanks to (|3.48p . Moreover by the 
following estimate of the chemical potential 



a 1 



2 7T V 2ttHo 

which is a consequence of p.22p . (|3.41[) and p.48|) . one obtains 



(3.53) 



[-|A + tt 2 £/(x) + 2 £ - 2 / 2 up (x) - fj%] f sup (x) > [ix" 1 ^ (2s - 1 + (ano)- 1 e 2 ) + 
an e- 6 (l - i (aftor 1 £ 2 - x) + 2e~ 2 / s 2 up (x) - M £ - 0(1)] /„„(*) > 
V^oe -2 (l-| (afl^e 2 - xj + 7r~ 3/2 exp | -erf} (l - (afi ) _1 e 2 -x) 2 j - 

2-1/2^-3/2 A + o(Q- 1/4 ))] f sup (x). (3.54) 



-(afl )- 1 e 2 -x) }> 



Now we estimate for any Si such that 1 — Si(aCl) 1 / 2 < 1 — (aflo) 1 e 2 

yfan Q e- 2 (l - (aftor 1 e 2 - x) + 7r~ 3/2 exp | - aft (l 

i (ar! )" 1/2 + ^ 3/2 - ^ 3/2 att (l - (afto)" 1 e 2 - xY > (1 - 5 2 )tt- 3 / 2 - CO 
where we have used the trivial inequality e~ x > 1 — x. On the other hand if x < 1 — (5i(aO) -1 / 2 



-1/2 
' 



v/a^£~ 2 (l - i (OA.)" 1 e 2 ~ <c) > h ( 



i-cn, 



-1/2 



Hence for any <5i such that 



1 1-2 -1 
— r < Sf < 1 =, 

2tt 3 V2 



(e.g., 5i = |) and fio sufficiently large, the r.h.s. of (|3.54|) is positive, i.e., / sup is a supersolution to (|3.47D 
and (|3.50[) follows. 

We now focus on the pointwise lower bound and prove that for any x £ such that 



|1 - x\ > 8 2 (afl)~ 1/2 , with log ^2 < S\ < 1, 



one has the lower bound 



5 2 » > (l-0(e 00 ))i7r- 3 / 2 v /oA^cxp{- a O(l-x) 2 }. 



Define, for any x £ 



i7T- 3 / 2 VoI^ 



exp |— ctfi (1 — x) 2 j — exp { — aQr,r] 2 } 
exp {— aSIS 2 } — exp {— afio?? 2 }] , 



, for |1 — a: I > S e 
otherwise, 



(3.55) 
(3.56) 

(3.57) 
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where S e := S(afl)- 1 / 2 and S = 0(1) is a suitable parameter. Clearly / sub vanishes at x = 1 — e 2 r] £ and, 
using (|3.53[) , one can show that / su b is a subsolution to p.40j) in A Vr : If |1 — x\ > 5 £ 

- ±A/ sub + n 2 Uf suh + 2e- 2 f sub 3 - ^ v / sub < 



^an + Tr-^ 2 ^n exp{-an(l-x) 2 }-^ + O(e- 2 Ve ) f suh < 



2- 1 / 2 7r- 3 / 2 v ^e- 4 [-1 + V2e~ 
if S 2 < log v2 and fig is large enough. On the other hand for |1 — x\ < S E 



-1/4 



/ S ub<0 (3.58) 



±A/ sub + n 2 Uf suh + 2e- 2 f suh 3 - ^v /sub < i anS 2 + 7r -3/2 v /^ e -* _ At gv + 0(£ -2j /sub < 



1 + Cfi 1/2 1 / sub < 0, (3.59) 



if 5 < 1 and ilo is large enough. Hence it suffices to take any S satisfying log s/2 < 5 2 < 1 to obtain that 
/ sub is a subsolution and g^ > / sub inside -A?;,,. However since for any x G ^ v /fj7, 

exp {§afi(l - xf - ±qA)?7?} < exp {-|aQ %(% - 1)} = 0(e°°), 



thanks to ((33H1) and (|3"12)l . then (j3"36) is proven. 

Collecting (QPgJ) , (|330)l and (|336) . we obtain (JS32J). 



□ 



Another crucial property of is its approximate symmetry w.r.t. x = 1: The potential [/ is symmetric 
w.r.t. a; = 1 up to small corrections, so one could expect that the same is true for g Ve . We then introduce 
another reference profile <? S ym which is defined as the unique positive minimizcr of the one-dimensional 
energy functional 16 ! 



fn/] := 2* / dx{i(/') 2 + ^ 2 2 (l-x) 2 / 2 + £ - 2 / 4 }, 

-e 2 n e 



f G ^([1 - e 2 % , 1 + e 2 n e ]) :/ = /*, 2tt / dx / 2 (x) = 1 



where / G ® sym with 



We also define 

E s ^:= inf i? sym [/]=i? sym [5 S ym]. 

/e@ aym 

We collect some useful properties of g S ym in the following 

Proposition 3.6 (Symmetric profile g sym ). 
The profile g sym satisfies the variational equation 

~ 2^sym "f" 2 a "^"(^ — 2 ') Ssym + 2e 5sym — l^ ^ Ssyrj 



(3.60) 

(3.61) 
(3.62) 



(3.63) 



with Neumann conditions at the boundary g 1 m (l±e 2 r] £ ) = 0. Moreover g sym is symmetric under inversion 
w.r.t. 1, i.e., 

g S ym{x) = # S ym(2 - x), (3.64) 

for any x G [1 — e 2 ?/ £ , 1 + £ *i /ias a unique maximum point at x = 1. 



3 Note that £ s y m coincides with £ aux except for the different integration domain. 
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Proof. Most of the results have already been proven in a slightly different setting in Prop. [5731 so we skip 
their discussion. The symmetry of g S ym can be obtained by exploiting the invariance of the variational 
equation under inversion or by a rearrangement argument. □ 

The relation between g Ve and g S ym is discussed in the next proposition but we first need to state 
an estimate of the gradient of g Ve . We actually state and prove two different estimates that we use for 
different purposes: (|3.65[) below yields the L°°-bound one should expect from the scales of the problem, 
whereas (|3.66[) yields the expected rate of exponential decay by giving a direct access to the ratio g' Ve g^ 4 - 

Lemma 3.1 (Estimate for gL)- 

If Q = floe~ i with fig > 0, then uniformly in x € A, h , 

<C(anf 2 e- cn ^ 2 , (3.65) 

|<4 (a;) | <Cnlr,le- 2 g Ve (x). (3.66) 

Proof. The proof of (|3.65[) is identical to that of |R1[ Proposition 2.2] (actually Step 3 of that proof). 
Starting from the L°°-estimate for g Vc we can use the variational equation to bound Ag., )e . Combining 
this with the Gagliardo-Nirenberg inequality [N| Theorem 1] we obtain the desired bound. 

To prove (|3. 66[) we can assume that x < x maK so that g' > 0, since the case x > x max is identical. By 
integrating the variational equation (|3.40[) from 1 — s 2 rj £ to x and using Neumann boundary conditions, 
we get 



<?;» = f dt Ut-'g'nM + - 2e- 2 gl(t) &U 2 (t)] g Vt (t)} > 



1 „r 

2 



- cn' 9rie ( X ) / dt (i - t.y > -CQ's^ig^ (x), 

where we have used the monotonicity of g Vr for x < x max and the sup estimate g Ve < ^e 2 fif^, which is 
the analogue of (|3.13[) . □ 

Proposition 3.7 (Estimate for g, h — <j sym ). 
If VI = Hqe ~ a with flo > large enough, one has 

\g„M ~ 9sym(x)\ < C^ 1 ^ (3.67) 

uniformly inxg A Ve ■ 

Proof. Several arguments of this proof can be taken over from a similar analysis in [Rll Section 2]. We 
thus focus on the main new features and skip some details for brevity 
ft is more convenient to work with rescaled variables: We introduce 

g„M = (any 1/4 g Vc [(an) 1/2 (l-x)) 
g sy m(y) = (an)~ 1/4 g sym i^aVt) 112 (f - a;) J 
that are minimizcrs of the functionals 

£%\f\ = 2*[ dy{||Si,/| 2 +(^y 2 + o((an)- 1 /2))/a + (anor 1/2 f A ) (l + {any 1 ' 2 y) , 



£ 



sym [/] = 2tt / dy {i \d y f\ 2 + \y 2 f + (an y 1/2 f 4 } 
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where the integration is on the domain 

The quantity 0((af2) -1 / 2 ) in the potential term of £|J comes from the Taylor expansion of the original 
potential. This imprecise but concise form is convenient because the precise expression is not relevant. 
Strictly speaking we should take into account the behavior of this term at infinity (since the domain 
A„ e extends to the whole real line in the limit we are considering) but any unwanted increasing term is 
compensated by the exponential fall-off of g Ve , inherited from that of g„ c . We do not elaborate more on 
this point as very similar subtleties have been dealt with in |R1[ Sections 2 and 4.1]. In the sequel we 
also replace A r)e by the whole real line. Again, thanks to the fall-off of the functions under consideration, 
this is harmless. 
We have 

H 9 y + + 2 ^ S 2 ym) Ssym - \ Sym 9- S yn„ (3.68) 

where we have set 

and A sym is a chemical potential fixed by the normalization, satisfying 

A sym = £ sym + £ f dy~g 4 sym . 

JR 

Because of (|3.68p . g S ym is an cigenstate of the Schrodingcr operator 

tfsym := + l y 2 + 2£g 2 ym . (3.69) 

As it is also positive (being the unique positive minimizer of (|3.60[1 ). it must be the ground state |LL1 
Corollary 11.9]. On the other hand, we have for g Vr 

{~h d l + W + 2 ^l) 9vs = *9v* + Vl* ((any 1/2 ) , (3.70) 

where 

X = Ei:+£ f dygl 

JR 

and the symbol O-u (d) denotes a function whose norm in the Banach space % is bounded by d. The 
remainder in (|3.70H has two origins (further details on a similar problem can be found in jRll Proof of 
Theorem 2.1]): 

• The remainder in the Taylor expansion of the potential. As already explained, the exponential 
decay of g„ e allows to control the errors produced by the higher order terms in the expansion. 

• The fact that the energy density in £|J is integrated against the measure (1 + (aQ)~ 1 / 2 y)dy. Again, 
most terms are dealt with using the fall-off of g Vc . Note that there is also a term of the form d y g Vc 
(coming from the replacement of the 2D Laplacian of a radial function by a double radial derivative) . 
This one is controlled using Lemma 13. II 

We now write 

h ■= g ne - g sy m- (3-71) 

It is a simple matter to prove that 

\\h\\ L2{s) <Ce (3.72) 

by noting that both g„ e and g S ym are close to g osc . Similar arguments have already been used so we omit 
the details. The rest of the proof consist in iteratively improving (13.721) . 
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Taking the difference of equations (|3.68[) and (|3.70p we see that 

(ffsym _ A sy m) h = (X _ ^ym) + ^ (-3^ _ + ^ ((c^)" 1 / 2 ) . (3.73) 

We have 

A - A sym = E% - E sym + 1 ( [ dy ~gl - [ dy g s 4 ym ) 
and easy energetic arguments combined with Schwarz inequality yield 

\\-\ s y™\<c( a ny 1/2 + ci\\h\\ L2m . 

Similarly 

\\9ve (Ssym - 9riJ || L 2( R ) < INL 2 (R) ■ 

Thus, multiplying (|3.73p by h and integrating, we obtain 

(h\ H°y™ - A sym \h) < c {any 1 ' 2 + a \\hfy m , (3.74) 

which simplifies to 

(h\ H sym - A sym \h) < Cmax((afi)~ 1/2 ,£ 3 ) 

by using (|3.72[) . We thus have, combining the above and (|3.72p and recalling that <? sym is the ground 
state of _ff sym associated with the eigenvalue A sym 



9Ve 



(1 + 0{£)) g sym + Ou^ (max [{ati)~ 1/2 , £ 3 ] ) (3.75) 



where "H sym is the Banach space obtained by taking the closure w.r.t. the norm induced by H sym (which 
in particular dominates the harmonic oscillator norm), i.e., 

H/I&sym : = </l- ffsym |/}> (3-76) 

and the remainder is orthogonal to g S ym- Since both <? sym and g r)e are also normalized one obtains from 
([3~75]) the identity 

1 = (1 + 0{£)f + O (max {any 1 , £ 6 ) 

which implies that, since £ <C 1, 

0{£) = O (max {any 1 ,£ e ) , (3.77) 
and therefore I < {any 1 , thus proving from (|3.75p that 

\\h\\ nsyn , < C {any 1/2 . (3.78) 

To obtain (j3~67)l it only remains 

• to recall that we are dealing with ID functions, thus the 'H sym -norm that controls the ii rl -norm 
also controls the L°°-norm; 

• to scale back the variables. 

□ 
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3.3 Energy Asymptotics and Consequences 

In this section we prove the energy asymptotics in the regime fi ~ e -4 and derive some bounds on reduced 
energies that will lead to the proof of the absence of vortices in the next section. As usual the asymptotics 
is derived from upper and lower bounds to the energy, the former being quite simple to obtain (see below) 
but the latter requiring more work. 

The starting point is an energy decoupling already used in Section 2.2 (in the context of GL and GP 
theories the technique originates in jLMj ). 

Proposition 3.8 (Reduction to A Ve ). 

Defining the function u £ H 1 ^^) by setting for x G A rie 

V> GP (x) =: ^,(s)u(x)exp{i|fi]0}, (3.79) 

the following holds 

E GP >E%+£ Ve [u], (3-80) 

with B given by (|3.1[) arid 

£,.[«]:=/ dx 5 2 E {i|Vu| 2 -B.( iU ,V U ) + e"X(i_| u |2) 2 |. (3 . 81) 

Moreover u satisfies the following equation on A Vs for some A S R 



V (s^Vu) - ig*B + (M 2 - 1) u = Xglu. (3.82) 



Proof. We first use the positivity of the energy density to restrict the integration in the GP functional to 
the annulus A rje - A standard energy decoupling for which we refer to |GRY[ Proposition 3.1] or |CPRY1[ 
Proposition 4.3] leads to (|3.80[) . The variational equation satisfied by g Ve is used, as well as the Neumann 
boundary conditions which ensure that there are no boundary terms in the decoupling. The derivation 
of the equation (|3.82j) uses the same ingredients, we refer to |CRY1 Lemma 4.3]. □ 

In order to conclude the proof of Theorem ll.41 it suffices to prove a lower bound to £ Ve [u] that matches 
the trivial bound which is obtained by taking u = 1, i.e., using g, h exp{i[£\|i9} as a trial function for £ GP . 
This function has of course to be extended to the whole space and normalized. Using (|3.42[) to bound 
g 7U outside A nE it is not difficult to estimate the remainders due to these operations and obtain 

E GP < £fj +0(e°°). (3.83) 

To prove the corresponding lower bound we have to analyze in some details the reduced functional 
£ Ve [u]. As is common in the context of GP theory (see, e.g., |AABI IPMT1 ICRYl ICPRYlj ]). a crucial step 
is an integration by parts of the angular momentum term. For technical reasons it will be convenient to 
first replace B with 

B(x) := O (r - r" 1 ) e fl (3.84) 

i.e., replace [^J with ft. We will also reduce the integration domain to A^: where the crucial L°°- 
cstimate Q3.43P holds true. These considerations lead to the definition of two potential functions that are 
'anti-derivatives' of g^B vanishing at the boundaries of A^^: 

F 1 {x):=nf At (t-t-^gKt), F 2 (x) := -Q / + ^ At (t - t~ l ) g^t). (3.85) 
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The use of two potential functions is a novelty of the present paper: When x is smaller than x max we 
have W ± Fi = <7jL B and will use F\ as potential function, whereas F 2 will be used for x > x max . In a 
sense this means that we use a potential function that is discontinuous at x max - 

This apparently strange strategy is actually very useful to deal with the boundary term that the 
integration by parts of the angular momentum term produces. With these definitions we shall have to 
estimate a boundary term on the circle of radius cc max , i.e., in the middle of the bulk, where g„ e is large. 
This is much easier than estimating the term located on a boundary of where the density is small, 

but such a term would be produced by considering a single potential function as in |AAB[ ICRYj . The 
necessity of this approach is due to the particular situation under consideration. The density profile g r)e 
behaves as a Gaussian and the long tails of such a function impose the use of a larger domain (which 
is apparent in the fact that A Ve is much thicker than the characteristic length of the Gaussian) than in 
former situations |AAB1 ICRY1 ICPRYlj where the profile was of TF type. This requirement rules out the 
trick we used in [CPRY1| . Sect. 4.3, to avoid having to estimate any boundary term at all. 

We now turn to the proofs of two lemmas about the potential functions we just introduced. We 
first estimate the difference between F\ and F2 at x max , in other words the jump of our 'discontinuous 
potential'. This quantity will appear in the boundary term produced by the integration by parts. The 
second lemma is a pointwisc estimate showing that F\ and F2 can be controlled by the density, which is 
required for the treatment of the bulk term. 

Lemma 3.2 (Estimate for Fi - F 2 ). 

For any x £ and in particular for x = x max 

\F 1 {x)~F 2 {x)\<Cn r i 3 £ / 2 . (3.86) 

Proof. We compute 

Fi(x) -F 2 {x)=Q df (t-t- 1 ) g 2 Ve (t) 

and note that a Taylor expansion yields 

t-\ =2(t-l)+0(e%) 
on the domain we are considering. We thus have 



- f 2 (x) = 2n dt{t- i)<4( t ) + o(n otfe ) = 



2n 



dt (t - 1) (<4(t) - g 2 ym (t)) + 0(fi o »fe), (3.87) 



where we have used the L 2 -normalization of g, lE and the symmetry of g S ym to obtain the first and second 
error terms respectively. Finally 



dt (t - I) (gl_(t) - g 2 sym (t)) 



< cn v 3 e /2 



by using the crucial estimate (|3.67p and the upper bound on g ne (|3.44[) . Recalling that % 3> 1, p.86p is 
proved. □ 



In the following statement we will denote by xs the characteristic function of a set S. 
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Lemma 3.3 (Pointwise estimate for F\ and F 2 ). 

If ft = £IqE ~ 4 and ilo > 0,q with Clo = 0(1) large enough but independent of e, then for any x £ A^jg 

X[i-e'^. Xa ,M\F 1 (x)\ +X[x^,i+e'^(x)\F 2 (x)\ < a- 1 (l + Cr>o 1/4 ) gl(x), (3.88) 
where the constant C > does not depend on x. 

Proof. We can assume without loss of generality that x < x max because in the opposite case the proof is 
essentially the same. For any x <G A^g, p.43[) yields 

{F^l^ ^n- 3 / 2 y^nr]ri(l + Cn^ 1/4 ) f dt{l-t 2 )c^p\-afl{l-tf\ < 

( poo 
l + CO(7 1/4 J / dz zexp{- a nz 2 } < 

±7r- 3 / 2 a-V 2 nY 2 (l + C% 1/4 )exp{-ari{l-x) 2 } < a' 1 {l + 0^^^ g 2 c (x). (3.89) 

□ 

The two above lemmas will be combined with some information about the function u. Again, wc 
distinguish the ingredients needed for estimating the boundary term from those needed for the bulk 
term. The latter are quite simple actually: The bulk term involves the vorticity measure 

/x := curl(ut, Vit) = \i curl (mVji* — ji*Vm) , (3.90) 

a quantity that is easily related to the kinetic energy density |Vu| 2 . Note that in our setting we do not 
need to apply tools involving a construction of vortex balls to obtain a useful estimate. 

Lemma 3.4 (Bound on the vorticity measure). 

We have 

\n\ < |Vu| 2 (3.91) 

pointwise in A Ve . 

Proof. A straightforward computation shows that 

\i = — 21m [d\u (c>2m)*] 
where di = d Xi for i = 1, 2 and Im stands for the imaginary part. On the other hand 

\dw (d 2 u)*\ <i(|9iw| 2 + |9 2 «| 2 ), 

which completes the proof. □ 

The estimate of the boundary term will rely on the fact that the degree of u around the circle of 
radius x max is in a certain sense controlled by the kinetic energy. Note that we do not know yet that u 
does not vanish on this circle, so (|3.93[) below is not truly a degree estimate. Note also that, as will be 
clear from the proof, the estimate holds because the circle we are interested in lies in a region where g Ve 
is large enough. This is the main motivation for the trick of choosing a potential function discontinuous 
along this particular circle. If the potential function was continuous there, we would have to estimate 
the degree on a circle were g„ e is small (namely on one of the boundaries of A^g). This would be much 
more difficult, if not impossible. 
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Lemma 3.5 (Circulation estimate). 

Let R be a radius satisfying 

r = i + e>(rr 1/2 ). 

We have, for some constant C > and any 5 > 0, 



da (iu, d T u) 



OBr, 



< c 



Proof. We use a smooth radial cut-off function \ with support in [R, R] , for some radius 

R = R-cQ~ 1/2 

with c > 0. Obviously one can impose 



(3.92) 
(3.93) 

(3.94) 



x(R) 


= 1, 


x(R) 


= o, 


\x\ 


< 1, 


|Vx| 





(3.95) 



We then use Stokes' formula to obtain 
da (iu, d T u) 



0B L 



/ da x(iw, d T u) = — / dx V ± x ' (* w : ^ w ) + / dx 

JdB R JK? JR 2 



where the integrals on R 2 are actually reduced to the support of x- If we choose c small enough in (|3.94l) . 
a combination of (|3.43p and p.92[) shows that 



sl > cn 1 ' 2 



in this region. Indeed, because of the assumptions p.92[) and (|3.94[) . the region of interest is located at 
a distance of at most 0(n -1 / 2 ) from 1, i.e., in the middle of the bulk (see, e.g., Remark 3.1). We thus 
estimate, using also (|3.91[) and p.95[) . 

/ da (iu, d T u) < Cfl 1/2 [ dx|u||Vn|+ / dx \n\ < 

JdB R JR<x<R JR<x<R 

cf d^glMv^ + cn- 1 ' 2 f dx 5 2 jV U | 2 < 

CS I dx 5 2 jVu| 2 + C(T 1 f dx g 2 h \u\ 2 + CQ- 1/2 j dx 5 2 jVu| 2 . (3.96) 

•^W "^W 

There only remains to recall that <7r? e M = \ip | is normalized in L 2 (M. 2 ) to complete the proof. 



□ 

We are now able to conclude the proof of the energy lower bound. A corollary of the proof below is 
the following useful estimate that will allow to complete the proof of the absence of vortices in the next 
subsection. 

Corollary 3.1 (Estimates for reduced energies). 

Denoting 

F Ve [u]:= I dxgl{\\Vu\ 2 +e- 2 g 2 u (l-\u\ 2 ) 2 }, 
one has for larger than some CIq = 0(1) 

\£ 7h {u}\ + F v Au}<cn 2 v 3 £ 

for some finite constant C < oo. 



(3.97) 
(3.98) 
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Proof of Theorem \1.4\ and Corollarv \3.1\ The starting point is the energy decoupling of Proposition 
and the trivial upper bound (|3.83[) . We now bound from below the reduced energy £ rk [u] . 

As announced we first replace B with B. The remainder produced by this operation is estimated as 
follows: 

/ dxj4(fi-|nj)ar 1 e,.(iu,Vu) <\Cp [ dx g\ \u\ 2 + iC/T 1 / dx. 9 2 jV W | 2 < 

iC/T 1 / dx 5 2 jV W | 2 + ±C/3, (3.99) 
JA ne 

where we have used the normalization of tp GF and (5 > is a constant to be fixed later in the proof. We 
now reduce the integration domain to ^4 v /^r at the price of a second remainder (for a lower bound we 
neglect the positive terms in the energy density) 



■Ave \-A~yve 



<\[3 I dx^|B| a |«| 



hr 1 



f dx g 2 Ve \\7u\ 2 < I dx <4 \Vu\ 2 + 0(s c 



Here we have used the exponential smallness of \ip GP \ = g Ve \u\ outside A^g. This requires 

T] e > |l0ge|, 

which we are free to decide. After these two steps we have 



e gp > E sv 

'/£ 



f dxr4(i| Vu | 2 -nB-(iu,Vu) + e- 2 gl(l-\u\ 2 ) 2 } 

J A .ACT 



+ / dx^lilV^ + ^^l-H 2 ) 2 }-^- 1 / gl\Vu\ 2 -\p~0{e° 

We now turn to the main part of the proof, namely a lower bound to 

£^[u]:= f dxgl{t\Vu\ 2 ~nB-(i Ul Vu) + e- 2 gl(l~\u\ 2 ) 2 }. 

As a first step we perform the integration by parts we have been alluding to. Since 

if 1 - £ 2 v / ?77 < x < x max , 



V ± F 1 = n g 2 B, 

V ± F 2 = ng 2 B, 



if a^max < x < 1 + e 2 Jrf e 



(3.100) 



(3.101) 



we have 



dx.g 2 e [\ \Vu\ 2 - 9.B ■ (iu, Vit)} 

/ dx \gl |V U | 2 + / 
JA./tt- Jl 



e 2 ^/rJ^<x<x Tt 



dx Fi{x)n 



x max <-x<l-\-e 2 yfrfe 



dx F 2 (x)fi 



[Fi(x max ) - F 2 (x max )] 



da {iu,d T u), (3.102) 



m r 
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where fx is the vorticity measure defined in p.90[) . We then combine Lemmas 13.31 and 13.41 to control the 
bulk terms produced by the integration by parts when f2o > Qo, the constant appearing in Lemma [ 

dx ^1 |Vu| 2 + J dx.Fi (a)// + f _dxF 2 {x)fi > 

For the boundary term we use Lemmas 13.21 and 13.51 to infer 



\Fi(x mllK ) - F 2 (x max )\ 



We thus obtain 



da (iu, d T u) 



da (iu, d T u) 



< 



Cn^S- 1 +Cn r) 3 £ /2 (o~ 1/2 +d) J dx.g2jVw| 2 . (3.103) 



J dx g 2 c |i |Vu| 2 - QB • (m,Vu)j 



> 



(± - Cn r, 3 J 2 S - a- 1 (l + 1/4 )) ) j dx g% c | Vu| 2 - C^ifJ 2 ^ 1 



We now recall that a > 2, which is obvious from its definition (11.53[) . It is thus sufficient to choose 

5 = OT - V 3/2 

with a sufficiently small constant C to deduce from the above that, for some C > 0, 



J dx^{i|V U | 2 -ttB-( iU ,V U )}>C^ 



>c I dx^j Vw | 2 ~o{n 2 4). 



Going back to p,102[) we have 

£^[u] > C [ dxgl \Vu\ 2 + [ dx 9 \ (1 - | U | 2 ) 2 - 0(nW E ). 
Combining this with (|3.83[) and (|3.100|) and taking 

= 0(1) 

large enough, we obtain for Qq > CIq 

E% + 0(e°°) > E GP > E% +C [ dx 5 2 e |V U | 2 + f dx ^ (l - \u\ 2 f - 0(e°°) 0(Q 2 oV 3 E ), 

JA r]c JA, le 6 

if again % ^> | loge|. For definiteness we fix 

r? e = |loge| 3/2 (3.104) 



in the sequel, a choice that indeed satisfies the requirements (|3.36|) . Neglecting the positive terms in the 
lower bound concludes the proof of Theorem II .41 whereas keeping them proves Corollary 13. II □ 
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3.4 Transition to the Giant Vortex State 

The proof of the absence of vortices requires one more ingredient, namely an estimate on the L°°-norm 
of the gradient of u. It is proved in much the same way than the corresponding results in |CRY1 ICPRYlj 
and we shall be brief on the details. 

Lemma 3.6 (Gradient estimate for u). 

For any c > 0, let A c be the domain defined as 

A c := jx : 1 - c| loge| 1/2 fr 1/2 < x < 1 + c| logel 1 / 2 ^- 1 / 2 } . (3.105) 
// the constant c > is chosen small enough in the above, there exists a C > such that 

IIVuIL-^) < Ce^-\ (3.106) 
Proof. The proof uses the equation (|3.82[) and the Gagliardo-Nirenberg inequality [Nj Theorem 1] 

HVuHi. < C (IHI^ ||Au||££ + \\u\\ LOB ) 
as in |CRY[ Lemma 5.1], which leads to the following bound (we omit the details) 

+ C ||^ GP ||^ ( p) \\ti\\ L -<p) (l A l 1/2 + II^V^J^^ + l|B||^ (P) ) (3.107) 

for any domain T> C A„ c . We will reduce to the domain (|3.105p and bound the terms on the right-hand 
side of the above equation. The most stringent requirement is the validity of a proper upper bound on 
, that requires a lower bound for g, h . Using p.43|) one can see that we have 

9r h > Ce^- 1 

on A c provided c is small enough. Note that with our choice (|3.104p . i.e., rj E = | loge| 3 / 2 , A c C A^fj^ 7 so 
one can use Q3.43P on A c . For the sup estimate on tp GF we use (|3.6p . whereas Lemma [3TT1 yields 

b^9vM L ^ {Ac} <Cfi 2 £ - 2 |log £ | 3 . 

On A c the vector potential B is easily seen to satisfy 

I|B|Il~(^) < Cn e- 2 \loge\. 

Finally, adapting the proof of |CRYi Equation (4.28)], we find 

|A| < |£,,[tt]| + Ce- 3 T, h [u} 1/2 < Cn vl /2 e- 3 , 

where the second inequality is a consequence of Corollary 13.11 and thus uses the assumption that Qq is 
large enough. Plugging these bounds in (|3.107p we obtain the result. □ 

We emphasize that on the domain A c defined in (|3.105p in the Lemma we have the estimate 

9 Ve > Ce^~\ (3.108) 

We are now able to conclude the proof of the absence of vortices. The idea that we use, namely the 
combination of bounds like those in Corollary [3T] and a gradient estimate, was first introduced in [BBH . 
We recall the definition of .Abuik: 

Auik = {x£l 2 : l-c\ logel 1 / 2 ^ 1 / 2 < x < 1 + c| logel 1 / 2 ^ 1 / 2 } , c < 
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Proof of Theorem ] 1.31 The argument is by contradiction. Let us assume that |1 — |tt(x)| | > | loge| a at 
some point x G -Atmik for some power a > 0. Then Q3.106P means that 



ii-Mi>i|io g£ r, 



on 6 x,C*e^— |loger a nAu 



ilk- 



This implies 



f dx%(l-| U | 2 ) 2 >C£ c2Q - 2 |log £ |- 

J A , £ 

This is a contradiction with p.98[) because of the choices 

1/2 



l0g£| 3/2 , 



c < 



Indeed, we would conclude that a negative power of e is controlled by a power of | log er| . We thus obtain 
|1 — M| < | loge|~ a in Amik for any power a > 0, and Theorem 11.31 is proved. □ 



The estimate of the degree of ip is a consequence of Lemma 13.51 and of the energy estimates of 
Corollary 13.11 

Proof of Theorem \1.5l Taking a radius R satisfying (|1.58l) we first note that the pointwise estimate in 
(|1.56[) implies that ip GP does not vanish on OBr, so that its degree is indeed well defined. We then 
compute 



2^deg(^ GP ,a^) 



dB r , 



GP 



GP 



GPI 



0B R u 



d(T 8 T ( -r—r 



= 2n[n\-i 

the second term on the second line being the degree of u. Then 

dcr (iu, d T u) 



dcr — d T 

dB R U 



< 



dBr, 



(3.109) 



(3.110) 



where we have used that |u| is bounded above and below by a constant on OBr. It remains to combine 
(|3.93p and (|3 .98(1 (where we stick to the choice (|3.104[) r\ e = \ loge| 3 / 2 for definiteness and optimize over 
S) and the result is proved. □ 



4 Rotational Symmetry Breaking 

We complete in this section the proof of Theorem 11.61 The strategy of the proof is identical to one 
followed in |CPRY11 Section 5], so we often omit technical details. 

A symmetric vortex minimizer is a function of the form f n {x) exp{in&} where f n is real and hgZ. 
For any given n there exists a minimizing function /„ with energy E n but we are here interested in a 
global minimizer, so we denote by n a minimizer of E n w.r.t. n£Z. The associated minimizing profile 
fn{x) is radial and positive far from the origin. Moreover it has a unique maximum at x max and is 
increasing for < x < x max and decreasing everywhere else. 

As in [CPRY11 Proof of Theorem 1.6] the result is proven by a direct inspection of the second variation 
of the GP functional around a local minimizer, which in the case of a symmetric vortex / n (x) exp{m?9} 



Critical Speeds in the GP Theory - CPRY - March 27, 2012 



45 



becomes 

Q„[H]:= / dx{i|(V- i A n )S| 2 + 7 ^|S| 2 +4 £ - 2 ^|S| 2 -^|S| 2 } + 

Jr 2 1 > 

2e~ 2 5R / dx/ 2 S 2 exp{-2mz?}, (4.1) 
Jr 2 

where 5 £ i/ 1 (IR 2 ) and \x n is the chemical potential associated with f n , i.e., \i n = E n + £~ 2 ||/„|||. 

Proof of Theorem \1.6i We assume that the symmetric vortex with degree n is a global minimizer of the 
GP functional and then show that this yields a contradiction since one can find a function 3 which 
makes the quadratic form Q n negative. 

However we start by certain useful properties of fn- Since we have assume that ip GF = f n (x) exp{m$}, 
the whole analysis contained in Section 13.11 applies to f n too and a simple inspection shows that it can 
be extended to any 51 > e~ 4 . In particular one can prove the estimates 

E GP <0(n), \\f n \\l<0(e 2 n), l/^^I^C^expj-fi 1 / 2 !!-^}, (4.2) 

for any x £ K 2 , which imply that f n is concentrated in the annulus 

An := |x e M 2 : 1 - VL- 1 ' 2 ] loge| 2 < x < 1 + ft 1/2 | loge| 2 ]} , (4.3) 

i- e -: H/nlU 2 ^) = 1 — The energy estimates (|3.41|) and (|3.83|) can as well be extended to angular 

velocities much larger than e -4 , implying 

0{e°°)>\ J dx(nx- 1 -r> a; ) 2 /I>(l-o(l))(n-r! + 0(r! 1/2 |log£| 2 )) 2 , 

which yields fi = 0(l + o(l)). 

Now we can introduce the test function which coincides with the one used in jCPRYll Proof of 
Theorem 1.6], i.e., 

S(x) := (A(x) + B(x))e^ n+d ^ + (A{x) - B{x))e l ^- d)i> (4.4) 

A(x) — ( xd+1 ^' for - x ~ Xmax ' B{x) — j™ 2 ^"' for < a; < x max , g . 

10, otherwise, l^lax/ni otherwise. 

A somewhat lengthy computation yields 

Q n [Z]=4ir f dx x 2d+2 f n f n {2(d + l)ti n + 2nn - (d + 2)n 2 x 2 - 2j(d+l){i 2 W{x)+ 
Jo 

poo 

jn 2 xW'(x) - 4(d + l)e- 2 f 2 } + 2irn 2 d 2 x™ llx \ Ax x' 1 f 2 . (4.6) 

Now using the facts 

» = n(l + o(l)), x max = 1 + o(l), W{x)=o(l), W'(x) = o(l), (4.7) 

for any x £ A n , as well as the monotonicity of f n for x < x max , we obtain 

Q n [E] < -47rfi 2 d(l - (1)) f °" X dx x 2d+2 f r J n + C'n 2 d 2 < 

Jo 

- 2nn 2 d [(1 - o(l))/ 2 (x max ) - Cd] < 0, (4.8) 
for e sufficiently small and any finite d > since f 2 (x mllx ) = 0(e 2 il) 3> 1. □ 
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